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Exact results for conformational statistics of compact polymers are derived from the two-flavour 
fully packed loop model on the square lattice. This loop model exhibits a two-dimensional manifold 
of critical fixed points each one characterised by an infinite set of geometrical scaling dimensions. We 
calculate these dimensions exactly by mapping the loop model to an interface model whose scaling 
limit is described by a Liouville field theory. The formulae for the central charge and the first few 
scaling dimensions are compared to numerical transfer matrix results and excellent agreement is 
found. Compact polymers are identified with a particular point in the phase diagram of the loop 
model, and the non-mean field value of the conformational exponent 7 = 117/112 is calculated for 
the first time. Interacting compact polymers are described by a line of fixed points along which 7 
varies continuously. 

PACS numbers: 05.50.+q, 11.25.Hf, 64.60.Ak, 64.60.Fr 



I. INTRODUCTION 

Lattice models of loops have emerged as an important 
paradigm in two-dimensional critical phenomena. They 
allow for a determination of the scaling properties of dif- 
ferent types of random walks which are used to model 
conformations of different phases of polymers [Q. For 
instance, the solution of the O(n) loop model has lead 
to exact results for conformational exponents of swollen 
and dense polymers as well as polymers at the theta 
point j^]. The theta point is the tricritical point which 
governs the transition between the swollen and the col- 
lapsed phase of polymers in solution |l| . Examples of con- 
formational exponents are 7, which describes the scaling 
of the number of polymer conformations with the num- 
ber of monomers Af, and v, for the scaling of the linear 
size of a polymer, as measured by the radius of gyration, 
with M . Here we calculate for the first time the exact 
value of 7 for polymers on the square lattice, in the com- 
pact phase. Compact polymers completely fill the lattice 
and are of direct relevance to statistical studies of protein 
folding |j§. 

Further motivation for studying loop models comes 
from the Fortuin-Kastclcyn construction which maps 
many discrete spin models (e.g., Q-state Potts) to ran- 
dom cluster models. Since cluster boundaries in two di- 
mensions form loops this naturally leads to a loop model 
representation. This random geometrical description of 
two-dimensional lattice models then provides a setting 
in which a general theory of their scaling limits can be 



sought. It is one of the goals of this paper to outline a 
specific proposal for such a theory in the form of an ef- 
fective field theory of fluctuating loops. This field theory 
is constructed following the Coulomb gas recipe || with 
some important new ingredients added It describes 
the fluctuations of a random surface for which the loops 
are contour lines. 

Scaling limits of many (but not all) two-dimensional 
lattice models are described by conformally invariant 
field theories H§). This observation has lead to exact 
results for critical exponents and other universal quan- 
tities, and to a classification of critical points based on 
their symmetry properties with respect to the group of 
conformal transformations. An obvious question which 
is often difficult to answer is: given a particular lattice 
model how does one construct the conformal field theory 
of its scaling limit? Loop models provide examples for 
which the scaling limit can be constructed in a physi- 
cally transparent way. This is accomplished by mapping 
a loop model to an interface model, where the loops are 
simply equal-height contours. An explicit coarse grain- 
ing procedure is then implemented for the height model, 
and it leads to a well known conformal field theory - the 
Liouville field theory. 

Interesting examples of loop models are also provided 
by one-dimensional quantum models, spin chains in par- 
ticular, where loops appear as world lines of the spin. 
This mapping of spins to loops has recently been used 
to formulate very efficient numerical schemes for simu- 
lating spin chains and ladders. These loop algorithms al- 
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low one to simulate much bigger system sizes and lower 
temperatures than by using more traditional algorithms 
with local updates 10 1. The loop representation of quan- 
tum spin chains also gives an illuminating stochastic- 
geometrical view of their quantum fluctuations 1 1 1 . For 
example, the spin-spin correlation function is related to 
the probability that two points on the space-time lat- 
tice belong to the same loop. This insight might lead 
to a practical theory of plateau transitions in the In- 
teger Quantum Hall Effect, i. e., one that would allow 
for a calculation of the correlation length exponent and 
other universal quantities which have been measured in 
experiments. Namely, the Chalker-Coddington network 
model , which is believed to be in the same universal- 
ity class as the plateau transitions, was recently mapped 
to an SU(n — > 0) quantum spin chain |l3| . It remains 
to be seen if this spin chain has a tractable loop-model 
representation. 

In the bigger picture, loop models are of interest as 
simple examples where the fundamental constituents are 
non-local, extended objects as opposed to point-like ob- 
jects such as particles and spins. Fluctuating geometries 
of this sort are used to model flux lines in superconduc- 
tors, domain walls in magnets, and crystalline interfaces, 
to name a few experimentally relevant systems. 

The extended nature of loops turns out to have pro- 
found consequences when one attempts to write down 
an effective continuum description of these models, say, 
following Landau's dictum of expanding the free energy 
(Euclidean action) in powers of the order parameter and 
its derivatives. Namely, terms which are geometrical in 
origin and non-perturbative in nature, and hence can- 
not be inferred from symmetry arguments alone, appear 
in the action. On the other hand, exactly because these 
geometrical terms are present the values of the effective 
coupling constants of the field theory are completely de- 
termined, a rather remarkable occurrence. 

Usually in an effective description provided by a field 
theory, coupling constants are phenomenological param- 
eters fixed by auxiliary information about observable 
quantities, such as the response functions or the related 
correlation functions. The Coulomb gas approach to two- 
dimensional critical phenomena is an example of an ef- 
fective theory wherein the electromagnetic coupling con- 
stant (i. e., the "magnitude of the unit charge") is de- 
termined from an exact solution of the model; typically 
it suffices to calculate the exact value of a single critical 
exponent. Our construction of an effective field theory 
of loop models closely parallels the Coulomb gas method 
with the important difference that the coupling constants 
are determined without recourse to any exact information 
about the model. For the model at hand no such informa- 
tion is available anyway, and moreover there are indica- 
tions that the model is not exactly solvable jl4j . On one 
level our theory can be viewed as a trick that allows one 
to calculate critical exponents in two-dimensional loop 



models without doing the "hard work" of exactly solving 
the model. On a deeper level it shows that lattice models 
of loops lead to continuum theories that are geometrical 
in nature, i.e., devoid of any couplings that depend on 
the microscopic details. 
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FIG. 1. Phase diagram of the two-flavour fully packed loop 
model on the square lattice. The loop model is critical for 
loop fugacities < nb,n g < 2. Particular points in the 
critical phase map to previously studied models: 6V - equal 
weighted six- vertex model [16], DL - dimer loop model [15], 
4C - four-colouring model [28]. The dashed line is the fully 
packed loop model studied numerically in Ref. [14]; the point 
CP along this line corresponds to the problem of compact 
polymers. Finally, the dotted line is the loop model for which 
an effective Liouville field theory was constructed in Ref. [7]. 

Here we study in detail the two-flavour fully packed 
loop (FPL 2 ) model on the square lattice. This is a sta- 
tistical model which describes two flavours of loops that 
occupy the bonds of the square lattice, subject to cer- 
tain close packing constraints to which we shall return 
shortly. The phase diagram of this model is described by 
two variables, rib and n g , which are the loop fugacities 
of the two flavours; see Fig. [I]. The phase diagram of the 
FPL 2 model has three important features that we wish 
to emphasize from the outset: 

i) For loop fugacities that fall into the region < n\> , n g < 
2 of the phase diagram the model is critical, i.e., it ex- 
hibits a power-law distribution of loop sizes. The novel 
feature is that every point in the critical region defines a 
different universality class characterised by an infinite set 
of geometrical critical exponents. All previously studied 
loop models (e.g., Q-st&te Potts, O(n) models) exhibit a 
line of fixed points. 

ii) The effective field theory of the FPL 2 model in the 
critical region describes a fluctuating two-dimensional in- 
terface in five dimensions, which is characterised by three 
elastic constants. We calculate these three couplings ex- 
actly as a function of the two loop fugacities. It is im- 
portant to note that all previously solved loop models are 
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characterised by a single elastic constant. 
Hi) From the field theory of the FPL 2 model we calculate 
for the first time exact results for the conformational ex- 
ponents of compact polymers on the square lattice. Fur- 
thermore, a particular line of fixed points in the phase 
diagram of the FPL 2 model can be identified with in- 
teracting compact polymers (rib = 0,n g < 2). We find 
that along this line the exponent 7 changes continuously, 
whilst v stays constant. 

The organisation of the paper is as follows. In Sec. |l| 
we review the scaling theory of compact polymers which 
provides our main motivation for introducing the two- 
flavour fully packed loop model on the square lattice in 



Sec. Ill . The rest of the paper is devoted to the study of 
this model using field theoretical techniques and numer- 
ical transfer matrix calculations. 

The FPL 2 model is mapped to an interface model in 



Sec. IV. For the interface model we construct the scaling 
limit in terms of a Liouville field theory, in Sec. |y|. In 
Sees. VI and VII we make use of the field theory to calcu- 



late the central charge and the infinite set of geometrical 
exponents associated with loops, in the critical region of 
the loop model. A short description of the non-critical 



region based on the field theory is given next in Sec. VIII. 

Following the field theoretical treatment of the FPL 2 
model, in Sees. IX and [X] we describe the construction 
of transfer matrices for different boundary conditions. 
They are used to determine the central charge, the first 
few geometrical exponents, and the residual entropy; the 
numerical results are in excellent agreement with the the- 
oretical predictions. Finally, in Sec. XI, we present some 
general observations regarding compact polymers and the 
Coulomb gas description of conformal field theories. We 
also comment on the dimer-loop model fl5|| and the three- 
state Potts antiferromagnet fl6| , in light of our solution 
of the fully packed loop model on the square lattice. The 
appendices are reserved for the calculation of scaling di- 
mensions of operators in the Liouville field theory and 
the enumeration of connectivities which are used for con- 
structing the transfer matrices. 



II. COMPACT POLYMERS 

Compact polymers, or Hamiltonian walks, are self- 
avoiding random walks that visit all the sites of the un- 
derlying lattice; see Fig. 0. They have been used as sim- 
ple models of polymer melts ]17|] and appear in statis- 
tical studies of protein folding ^,|). Unlike dilute and 



dense polymers whose scaling properties were calculated 
exactly from the 0(n) loop model JTsI , compact poly- 
mers defied a similar treatment until recently. Numerical 
transfer matrix calculations [ [l~g| ] , a Bethe-ansatz solution 
p0| , and a Coulomb gas theory |2l|] of the fully packed 
loop model on the honeycomb lattice, all conclude that 
compact polymers define a new universality class of crit- 
ical behaviour. Here we study compact polymers on the 
square lattice. We calculate exact scaling exponents and 
find them to be distinct from the honeycomb case. This 
was first reported in Ref. |L4| on the basis of numerical 
transfer matrix results. 
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FIG. 2. Compact polymer on the square lattice; x and y 
are the positions of the chain ends. 

The lattice dependence of critical properties distin- 
guishes the compact polymer problem from its dilute and 
dense counterparts in a crucial way. It places them into 
the class of geometrically frustrated critical systems 1 . 
A physically relevant measure of frustration for com- 
pact polymers is the number of contacts per monomer. 
Contacts are realised by monomer pairs where the two 
monomers are nearest neighbors on the lattice but are 
not adjacent along the polymer chain. In lattice models 
of proteins hydrophobic interactions among the amino 
acids occur at contacts {|,[| . For the square model stud- 
ied here the number of contacts per monomer is two, 
whilst on the honeycomb lattice it is one. 

In order to study the scaling properties of compact 
polymers we focus our attention on the two most widely 
studied conformational exponents v and 7. If R is the 
radius of gyration of the polymer then 



(2.1) 



1 Another example is the antiferromagnetic three-state Potts 
model which has a zero-temperature critical point on the 
square and the Kagome jz| lattices characterised by dif- 
ferent critical exponents. 
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where Af is the number of monomers. Since compact 
polymers visit all the sites of a lattice, they are space- 
filling and we conclude that v = 1/2. This simple result 
will serve as an important check on our field theoretical 
calculations where it will be recovered. 

In order to define the conformational exponent 7 we in- 
troduce C(Af), the number of compact polymers (Hamil- 
tonian walks) on a square lattice with Af sites. Since 
a compact polymer fills the lattice, boundary conditions 
(free, periodic, etc.) play an important role. Following 
Saleur and Duplantier |23| ], we define 7 in a way that is 
insensitive to the boundaries. Namely, if we introduce the 
quantity C a {Af), the number of compact-polymer rings, 
then we can expect 



7 = 117/112 = 1.0446- 



(2.5) 
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(2.2) 



where 7 does not depend on the choice of boundary condi- 
tions. Therefore, in order to calculate 7 we need to solve 
the hard combinatorial problem of counting the num- 
ber of open and closed compact polymers on the square 
lattice. Following de Gennes we do this by mapping the 
counting problem to the calculation of a correlation func- 
tion in a particular statistical model at the critical point. 

Consider the quantity Z(x, y;Af), the number of com- 
pact polymer conformations that start at the vertex x of 
the \fW x square lattice, and end at y (see Fig. ||); 
we consider the limit 1 <C |x — y| <C VAf, where x and 
y are chosen far from the boundaries of the lattice. For 
this quantity we can write down the scaling form |23|] : 



Z(x,y;A0 =C o (A0|x-y|- 2a;i / 



|x-y| 
TV 1 / 2 



(2.3) 



where f(u) is a scaling function with the property f(u) — ► 
const, as u — ► 0, and x\ is a geometrical exponent related 
to 7. Integrating Z(x,y;Af) over all end-points y and 
comparing the result to Eq. (2.2), the scaling relation 



7 = 1 - Xi 



(2.4) 



follows. 

To calculate the geometrical exponent x\ we introduce 
the two-flavour fully packed loop model on the square 
lattice. The fact that we need two loop flavours follows 
from the simple observation that the bonds not covered 
by the compact polymer also form loops whose number is 
unconstrained. For the loop model we then construct an 
effective field theory in which Z(x, y; Af) becomes a two- 
point correlation function. The asymptotics of this func- 
tion can be calculated exactly and we find x\ = —5/112, 
from which 



follows. This is to be compared to the mean- field theory 
value 7mf = 1 pij , which is also the result obtained for 
compact polymers on the honeycomb lattice p0[ . 

The conformational exponent 7 was measured di- 
rectly from enumerations of conformations of chains with 
lengths up to 30 in Ref. Qj, and the value 7 = 1.01(5) was 
reported. More recently, from a numerical transfer ma- 
trix study of the fully packed loop model on the square 
lattice the geometrical exponent x\ = —0.0444(1) was 
determined jlS, in excellent agreement with the exact 
result. 

Another quantity of interest is the connective constant 
k which determines the leading, exponential with system 
size scaling of the number of compact polymers |25| 



C{AT) ~ K M Kf d ~ iy 'AT- 1 



(2.6) 



Here k s is the surface connective constant; it appears due 
to the space-filling nature of compact polymers. Both the 
value k — 1.475(15) found in Ref. and the estimate 
k ~ 1.472 obtained from transfer matrix calculations sim- 
ilar to ours [ p6[ , seem in favour of the mean-field result 
Kmf = I = 1.4715 • • • [||. 2 In Sec. [XD| we report the 
very accurate numerical value 



k = 1.472801(10), 



(2.7) 



which shows that the connective constant for compact 
polymers also deviates slightly from the mean-field re- 
sult. 

For the remainder of the paper we elaborate on the 
calculation of 7 for compact polymers, in the process 
unveiling an extremely rich phase diagram of the asso- 
ciated loop model. As remarked earlier, it contains a 
two-dimensional region of fixed points, which we char- 
acterise in detail by calculating the central charge and 
the geometrical exponents associated with loops for each 
point on the critical manifold. 



III. TWO-FLAVOUR LOOP MODEL 

The two-flavour fully packed loop model on the square 
lattice was introduced in Ref. |^8| as the loop represen- 
tation of the four-colouring model pjJ]. It is the nat- 
ural generalisation of the fully packed loop model on 
the honeycomb lattice, which is the loop representa- 
tion of the three-colouring model pi] ]. In general, a q- 
colouring model on a g-fold coordinated lattice is given 



2 Very recently the field theory of Ref. |2^] has been improved 
p7) yielding, however, unchanged values for 7mf and k M f- 
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by edge colourings of the lattice with q different colours; 
an edge colouring of a graph is one where no two bonds 
that share a common vertex are coloured equally. The 
colouring model is mapped to a loop model by choosing 
[q/2] colour-pairs; each pair defines strings of alternating 
colour that necessarily form loops (unless they terminate 
at the boundary). In this way we end up with a loop 
model with [q/2] flavours of loops. 

To define the FPL 2 model we first specify the al- 
lowed loop configurations Q. In Q every bond of the 
square lattice belongs to one and only one loop of either 
flavour, and loops of the same flavour are not allowed 
to cross. Representing the two flavours by solid (black) 
and hatched (grey) line segments respectively this fully 
packing constraint allows each vertex of the square lat- 
tice to have one of the six appearances depicted in Fig. ^. 
Each loop is assigned a fugacity depending on its flavour: 
?ib for black loops and n g for grey loops. The partition 
function of the FPL 2 model is then 



Z 



(3.1) 



The fully packed loop model of Batchelor et al. is 
obtained by setting the loop fugacity of the grey loops 
to unity. In the limit — ► we recover the compact 
polymer problem. 



n a 



r 



J. 



FIG. 3. The six vertex configurations of the FPL 2 model 
that are allowed by the fully packing constraint. Black and 
grey loop segments are shown here as solid and hatched lines 
respectively. Each vertex is adjacent to four edges, here shown 
as dots, that are referred to as "dangling" if they are not con- 
nected to an edge of a neighbouring vertex. Note that the two 
rightmost vertices explicitly permit the two flavours to cross. 

If we define a restricted partition function of the FPL 2 
model, to which only configurations with a single black 
loop segment propagating between points x and y con- 
tribute, then Z(x,y;Af) in Eq. ( |2.3| ) is obtained in the 
limit rib —* 0, n g — > 1. The first limit discards all con- 
figurations with black loops present, leaving only the 
black Hamiltonian walk (compact polymer) between x 
and y, whilst the second ensures that all walks are 
weighted equally. We could also consider weighting dif- 
ferent Hamiltonian walks differently by setting n g ^ 1. 
This situation can be interpreted as describing interact- 
ing compact polymers, and, as will be shown later, it 
leads to a continuously varying exponent 7. A similar 
property of interacting oriented polymers in the swollen 
phase was suggested by Cardy from a field theoretical 



calculation |30|] . Recent numerical studies of the inter- 
acting oriented self-avoiding walk by Trovato and Seno 
plj, though, seem to be at odds with Cardy's predic- 
tion of an exponent 7 that varies continuously with the 
interaction strength. 

Some idea of the phase diagram of the FPL 2 model 
as a function of and n g can be gotten by examining 
the extreme limits of the loop fugacities. Namely, for 
n b, n g — > 00 all loops have the minimum length of four, 
i.e., they each surround a single plaquette of the square 
lattice. There are no large loops in the system and the 
model is non-critical, or in other words, the average loop 
length is finite. On the other hand, in the critical phase 
of the loop model, which is the subject of this paper, in 
a typical configuration one finds loops of all sizes char- 
acterised by a power-law distribution. This leads to an 
average loop length which diverges with the system size. 
Such is the case in the other extreme limit of loop fugac- 
ities, rib, n g — > 0, when the loops cover the whole lattice. 

Other previously studied models that are particular 
points in the phase diagram of the FPL 2 model are 
the four-colouring model, the dimer loop model, and 
the equal weighted six-vertex model; see Fig. |l|. For 
(rib,n g ) = (2,2) the loop fugacity of each loop can be 
evenly (1+1) distributed among the two ways of colour- 
ing the bonds occupied by the loop with two colours in 
an alteranting fashion: ABAB . . . for black loops and 
CD CD . . . for grey loops. This is then the symmetric 
four-colouring model (A,B,C, and D are the colours) 
studied by Baxter p2fl . In the dimer loop model black 
and white dimers are placed on the square lattice so 
that every vertex is covered by one of each |l5f] . If we 
identify the dimer covered bonds with the black loops 
then this model is mapped to the (rib,%) = (2, 1) FPL 2 
model. And finally (nb,n g ) — (1,1) constitutes the 
equal- weighted six- vertex model Q , the allowed vertices 
being those of Fig. |[ 



IV. HEIGHT REPRESENTATION 

The critical phase of the FPL 2 model can be described 
in terms of an effective field theory, following the general 
procedure discussed in Ref. ||. The idea is to think of 
loops as contours of a scalar field, which we refer to as 
the height. Depending on the loop model in question the 
height can have one or more components. If the number 
of components is D±_ then the effective field theory of 
the loop model describes a fluctuating two-dimensional 
interface in D± + 2 dimensions. 

To introduce the heights we first map the loop model 
to an oriented loop model, as shown in Fig. ^. The 
orientation of every loop is chosen randomly and inde- 
pendently. Every non-oriented loop configuration is thus 
transformed into an oriented one (G')', the number of 
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oriented configurations that correspond to the same non- 
oriented loop configuration is simply 2 Nb+Ns . 
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FIG. 4. Mapping of the FPL 2 model to an interface model, 
(a) — > (b): Transform the loop configuration into an oriented 
loop configuration by choosing the orientation of each loop 
independently and randomly, (b) — > (c): Every bond in the 
oriented loop configuration is in one of four states, depending 
on its flavour and direction; these four states are represented 
by three-vectors A,B,C, and D. (c) — > (d): The micro- 
scopic height z of the interface model changes from plaquette 
to neighbouring plaquette by A, B, C, or D depending on the 
state of the bond between the two plaquettes. The change in 
z is positive going clockwise around even vertices and coun- 
terclockwise around odd ones. 

Next, for each loop we redistribute its weight (fugac- 
ity), rib or n g depending on whether it is black or grey, be- 
tween the two possible orientations. For the black loops 
we do this by assigning to, say, the clockwise orienta- 
tion the phase factor exp(i7reb), and the opposite phase, 
exp(— i7reb), to a counter-clockwise oriented black loop. 
Similarly for grey loops the clockwise oriented ones are 
assigned a weight exp(i7re g ) whilst the counter-clockwise 
loops are weighted with exp(— i7re g ). The loop fugacities 
are related to the newly introduced parameters eb and 
e g by 



n h 



2 cos(7reb) 
2 cos(7re g ) 



(4.1) 



since the partition function of the original (non-oriented) 
model, as given by Eq. (3.1), must be recovered by inde- 
pendently summing over the two possible orientations for 
each loop. Note that for < nb,n g < 2 the parameters 
eb and e g are real, whilst for n\,,n g > 2 they are purely 
imaginary. As discussed in more detail in Sec. VIII this 
is the crucial property that leads to a critical state of the 
loop model in the former and a non-critical one in the 
latter case. 

Now that the loops are oriented we can interpret them 
as contours of a height field; the orientation is necessary 



as it determines the direction of increasing height. The 
systematic construction of the microscopic heights sets 
out from the observation that every bond of the square 
lattice is in one of four possible states: it can be coloured 
black or grey, and oriented from an even to an odd site, 
or from odd to even. "Even" and "odd" refer here to the 
two sublattices of the bipartite square lattice; every even 
site is surrounded by four nearest neighbouring odd sites, 
and vice versa. 

The four possible bond-states are represented by four 
vectors - which are the colours in the four-colouring rep- 
resentation - A,B,C and D; see Fig. |]c. The micro- 
scopic heights {z} are defined on the dual lattice and 
the change in height when going from one plaquette cen- 
tre to the next is given by A,B, C or D, depending on 
the state of the bond which is crossed; Fig. |Jd. For the 
height to be uniquely defined the four vectors must sat- 
isfy the constraint A + B + C + D = 0. This means that 
the microscopic heights live in a three-dimensional vector 
space, which we take to be Z 3 . In other words the ori- 
ented FPL 2 model maps to a model of a two-dimensional 
interface in five spatial dimensions. 

By reasons of symmetry the four vectors are chosen so 
as to point from the centre to the vertices of a regular 
tetrahedron. With a suitable choice of coordinates they 
are represented by three-vectors: 



A = (-1,4-1,4-1) 

B = (+1,+1,-1) 

C = (-1,-1,-1) 
D = (+1,-1, +1) . 



(4.2) 



This is the same normalisation as the one used in 
Ref. ||. 

Mapping the loop model to on oriented loop model also 
allows for a local redistribution of the loop weights. This 
is important since it leads to a local field theory for the 
heights. As we will find out shortly, though local, this 
field theory is somewhat unconventional due to the non- 
local, extended nature of the fundamental microscopic 
objects it purports to describe. 

To redistribute the phase factors associated with ori- 
ented loops we assign a phase exp(— «7reb/4) to a vertex 
of the square lattice if a black loop makes a left turn at 
that vertex, the opposite phase exp(+i7reb/4) if it makes 
a right turn, and the weight 1 if it continues straight. 
The total vertex weight A(x) is a product of the phase 
factor originating from the black loop and an equivalent 
one from the grey loop passing through the same v erte x 
x. The partition function of the FPL 2 model, Eq. ([O]), 
can now be rewritten as a sum over oriented loop config- 
urations (i.e., colouring configurations) 



EIIa(x) 

gi x 



(4.3) 
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Once the height at a single point is fixed Q' is in a one-to- 
one correspondence with the configurations of the micro- 
scopic heights, and the summand in the above equation is 
the appropriate weight. In the critical phase of the FPL 2 
model the interface described by Eq. ( |4.3| ) is rough, and 
the field theory is constructed so as to correctly repro- 
duce its long-wavelength fluctuations. 



A. Spectrum of electromagnetic charges 

The mapping from oriented loop configurations, which 
are equivalent to edge colourings, to microscopic height 
configurations is one to many. In particular, two height 
configurations corresponding to the same edge colour- 
ing can have their heights shifted with respect to each 
other by a global shift meK. The set TZ forms a three- 
dimensional Bravais lattice, i.e., it is closed under integral 
linear combinations, and its elements are the magnetic 
charges in the Coulomb gas representation of the FPL 2 
model. The lattice reciprocal to the lattice of magnetic 
charges, TZ*, defines the electric charges e £ TZ* , with the 
property e • m = 2itm, m S Z. 

The construction of the lattice TZ for the FPL 2 model 
follows the usual prescription for height models, and has 
been carried out in detail in Ref. p£|. For the sake of 
completeness we outline this construction below. 

It is convenient to first identify the flat states (also re- 
ferred to as the ideal states), i.e., those colouring states 
which minimise the variance of the microscopic height z. 
From the height mapping described above it follows that 
these states have all of their plaquettes coloured with two 
colours only; an example is shown in Fig. ||c. This leads 
to a colouring state that is periodic, with the same 2x2 
colouring pattern repeated throughout the lattice. There 
are twenty four flat/ideal states for the colouring repre- 
sentation of the FPL 2 model, corresponding to the num- 
ber of permutations of four different colours. Namely, an 
ideal state is completely specified by listing the colours 
of the bonds around a single site (say the origin), start- 
ing from the left horizontal bond and proceeding clock- 
wise. To each fiat state we assign a coarse grained height 
h = (z), which is the average microscopic height over a 
2x2 unit cell of the colouring. 

The flat states form a three dimensional graph, which 
we refer to as the ideal state graph, X. Namely, start- 
ing from any ideal state four other ideal states can be 
reached by exchanging a pair of colours that form a pla- 
quette. For example, by exchanging the colours A and 
B in Fig. all the ABAB plaquettes are turned into 
BABA plaquettes to give a new ideal state. Under these 
plaquette flips only the microscopic heights at the centres 
of the affected plaquettes are changed. In this way the 
ideal states form a four-fold coordinated graph in height 
space, where each vertex is indexed by a colour permuta- 
tion, and its position in IR 3 is given by the coarse grained 



height h. Bonds are associated with transpositions of 
two colours; they lie along the direction defined by the 
difference of the two colour vectors, and have a length of 
\/2/2 if the normalisation in Eq. (4.2) is chosen. 

The ideal state graph is a tiling of IR 3 with truncated 
octahedra; this regular polyhedron is better known as 
the Wigner-Seitz cell of a body-centred cubic (bcc) lat- 
tice (see Fig. ||) . A single truncated octahedron in T has 
twenty four vertices corresponding to the twenty four dif- 
ferent ideal states. The set of vertices in X representing 
the same ideal state form the repeat lattice TZ, which is 
face-centred cubic (fee) with a conventional cubic cell of 
side 4. 




FIG. 5. The ideal state graph of the FPL 2 model in the 
four-colouring representation. 

To obtain the continuum description of the FPL 2 
model we coarse grain the microscopic height over do- 
mains of ideal states. This gives rise to the coarse grained 
height h which we can consider to be compactified on 
IR 3 /7?.. The phase space of the height is not simply con- 
nected, thus allowing for topological defects (vortices) 
with topological charges that take their values in TZ [[34j . 
These defects are associated with magnetic charges in the 
Coulomb gas representation of the FPL 2 model. Electric 
charges on the other hand are associated with vertex op- 
erators exp(ie • h). If we take the height to live in IR 3 /7?. 
then vertex operators are well defined only for values of 
the electric charge e 6 TZ* . TZ* is the lattice reciprocal 
to the lattice of magnetic charges TZ, and it is a body- 
centred cubic (bcc) lattice with a conventional cubic cell 
of side 7r. 



V. CONSTRUCTION OF THE FIELD THEORY 

An effective field theory of the FPL 2 model should de- 
scribe large scale properties of loops. The kind of ques- 
tions we expect it to answer are ones that do not refer 
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to the microscopic details of the lattice model. For ex- 
ample, from the effective field theory we will calculate 
the asymptotics of the probability that two points lie on 
the same loop, when the separation between the points 
is large compared to the lattice spacing. From this and 
related quantities the conformational exponents of com- 
pact polymers can be extracted. 

The field theory of the FPL 2 model is defined by the 
Euclidean action for the coarse grained height h. Con- 
sider a typical configuration of the oriented FPL 2 model 
which is equivalent to the colouring model. It consists 
of domains of ideal states. To each ideal state domain 
we assign a coarse grained height, defined earlier as the 
average microscopic height over the domain. In the con- 
tinuum limit we assume that this height is a smoothly 
varying function of the basal plane coordinates (x l ,x 2 ). 
The partition function that takes into account only the 
large scale fluctuations of the height can be written as a 
functional integral, 

Z > = Jvh exp(-5[h]), (5.1) 

where S is the Euclidean action of a Liouville field the- 
ory with imaginary couplings Q|. The Liouville action 
contains three terms, 

S = S E + S B + S L . (5.2) 

Each one has a concrete geometrical interpretation in the 
FPL 2 model, which we describe next. 

A. Elastic term 

The first term in the effective action for the FPL 2 
model describes the elastic fluctuations of the inter- 
face. It gives less weight to configurations that deviate 
from the flat states, by penalising finite gradients of the 
height. This term is entropic in origin. Namely, in order 
to change the colour of a particular bond in the four- 
colouring representation of the loop model, say C — > B, 
all the C's and B's have to be interchanged along the 
CB loop which contains the chosen bond. This transfor- 
mation we call a loop flip; see Fig. |(| The ideal states 
maximise the number of loops of alternating colour and 
consequently they have the largest entropy of loop flips. 

d q± one 
oot^ oU[ 

n n r m m r 

FIG. 6. A loop flip changes one oriented loop configuration 
into another. Here the bond states C and B are exchanged 
along a single BC plaquette (cfr. Fig. ^c). 



In its most general form the elastic term in the effective 
action can be written as a gradient expansion, 

5 E = ~ J ePxK^ fi a i h a d j h p , (5.3) 

where higher powers of the height gradients and higher 
derivatives of the height are less relevant at large scales. 
The stiffness tensor K^a nominally has 36 components; 
the indices i,j — 1,2 are for the basal plane coordi- 
nates, whilst a, /3 — 1,2,3 label the three components of 
the height. Summation over repeated indices is assumed 
throughout. 

The number of independent non-zero components of 
the stiffness tensor (i.e., elastic constants) is actually only 
three, once all the symmetries of the FPL 2 model are 
taken into account. The relevant symmetry transforma- 
tions, that is the ones that become the symmetries of the 
effective action, are the ones that leave the weights of 
oriented loop configurations unchanged. First, there are 
the lattice symmetries, translations and rotations, which 
cut the number of independent elastic constants down to 
six. The terms that are allowed in Se are scalars under 
rotations in the basal plane {(x 1 , x 2 )}, and they are nec- 
essarily of the form dh a ■ dh 13 ', where d = (81,82) is the 
usual gradient. Second, the FPL 2 model possesses colour 
symmetries, 

A <-* B : eb <-> — e\> and Z\ «-> 2:3 (5-4) 

and 

C <-> D : e g <-» — e g and Z\ <-> — Z3 , (5-5) 

which interchange the colours and at the same time trans- 
form the microscopic heights. Taking into account the 
colour symmetries the elastic contribution to the action 
takes on the form: 

Se = \J d 2 x {K lx [{dh l f + (dh 3 ) 2 ) + (5.6) 
2K lz (dh l ■ dh 3 ) + K 22 (dh 2 ) 2 } . 

Furthermore, by introducing a change of coordinates in 
height space, 

H 1 = -{h 1 -ti i ), H 2 = h 2 , H 3 = ^{h l + h 3 ) (5.7) 

S'e becomes diagonal, 

S E = ~ J d 2 x g a (dH a ) 2 . (5.8) 

The three coupling constants g a (a = 1, 2, 3) are linearly 
related to the three elastic constants, 
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g 1 = 2(K 11 -K 13 ), g 2 = K 22 , g 3 = 2(K U + K 13 ). 

(5.9) 

The appearance of three elastic constants is rather in- 
triguing from the viewpoint of loop models that have 
been solved previously. The Q-state Potts, the O(n), 
and the honeycomb FPL models are all characterised by 
a single coupling constant, which has been determined 
case by case from their exact solutions. Below we will 
show that all three couplings in Eq. ( |5.8| ) can be calcu- 
lated exactly from the loop ansatz introduced in Ref. |7j . 3 
The ansatz states that the operator which enforces the 
complex weights assigned to oriented loops is marginal 
in the renormalisation group sense. This property of the 
field theory is intimately related to the random geom- 
etry of loops; we elaborate on this important point in 
Section |VC2|. 



to be non-zero there must be at least a single winding 
loop present. If this winding loop is black, then the height 
difference is A or B depending on its orientation; simi- 
larly if the loop is grey the height difference is C or D. 
Furthermore if the background charge is chosen so as to 
satisfy 



eo • A = 7reb eo • B = — net, 
eo • C = 7re g eo • D = — ire^ 



(5.11) 



then the winding loops will be assigned their proper 
weights. This is again seen by summing over the two 
possible orientations of the winding loop. In the normal- 
isation chosen for the colour vectors, Eq. (4.2), the unique 
solution of the system of linear equations in Eq. ( |5.11 ) is 



eo 



-(e g + e b ,0,e g - e b ). 



(5.12) 



B. Boundary term 



The mapping of the loop model to an oriented loop 



model with local complex weights A(x) (Eq. (4.3)) fails 



for loops that experience the boundary. For example, if 
we define the FPL 2 model on a cylinder then loops that 
wind around the cylinder will not be weighted properly. 
The winding loop has an equal number of left and right 
turns and hence it will be assigned a weight one. Sum- 
ming over the two orientations gives a weight two, and 
not the correct rib or n g , depending on the flavour. To 
correctly weight these loops one introduces a boundary 
term into the effective action, 



S R = 



^ J d 2 x (e • h)K ; 



(5.10) 



1Z is the scalar curvature and eo is the background elec- 
tric charge, which is to be determined. Since we are only 
concerned with the situation where the lattice on which 
the FPL 2 model is defined is flat, the scalar curvature 
vanishes everywhere except at the boundary. 

To determine eo we consider the FPL 2 model on 
the cylinder. The scalar curvature of the cylinder is 
proportional to the difference of two delta functions 
situated at the two far ends of the cylinder: 1Z = 
4tt [S(+oo) — S(— oo)]. Therefore Sb has the effect of 



placing vertex operators exp(±ieo • h) at 



= ±oo; 



here x is the coordinate along the length of the cylin- 
der. These vertex operators assign an additional weight 
exp(ie • (h(+oo) — h(— oo)) to oriented loop configura- 
tions on the cylinder. Now, in order for h(+oo) — h(— oo) 



This calculation of the vector background charge gener- 
alises the scalar case studied previously . 



C. Liouville potential 

The elastic term and the boundary term make up the 
usual Coulomb gas approach to two-dimensional critical 
phenomena. Recently we have argued that this descrip- 
tion is incomplete and that an extra term Sl must be 
added to the effective action. To see this consider a large 
loop in the bulk, one that does not experience the bound- 
ary. Without the extra term this loop would be weighted 
exclusively by the bulk term Se- There are two problems 
with this: Se is real whilst an oriented loop should be 
weighted by a complex phase, and, Se does not distin- 
guish between the two orientations of a loop which are 
assigned different weights. We conclude that an extra 
bulk term is necessary! 

The most general form of a bulk term is 



Si = d 



;[h(x)] , 



(5.13) 



wher e e xp(— w[h(x)j) is the scaling limit of A(x) in 
Eq. ( |4.3| ). In this sense Si, is energetic in origin, as op- 
posed to Se, which we argued in Sec. VA accounts for 
the entropy of edge colourings. 

Microscopically, the vertex weight A can be written in 
terms of the colours of the bonds around the particular 
vertex as A = exp(— w) where 



3 The coupling constant g for all the loop models known to 
date can be calculated using this method, therefore dispensing 
with the need for an exact solution. 
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w(B, C, A, D) 
u>(B,D,A,C) 
u>(A,B,C,D) 

iu(B,A,C,D) 

t«(A,B,D,C) 

w(B, A,D,C) 



= 0, 
= 0, 

= = Fi^(e g + e b ), 
• 7r / 

= Ti^(e g - e b ), 
= Ti^(eb - e g ), 
= Ti^l-eb - e g ) 



(5.14) 



the top sign is for even vertices whilst the bottom sign 
applies to odd vertices of the square lattice. Here we 
adopt the notation (ci, 02, 03, 0-4) for the ordering of the 
colours around a vertex by listing the colours clockwise 
from the leftmost bond. The operator w is completely 
specified by the values it takes on the six edge colour- 
ings listed above since it does not change under cyclic 
permutations of its arguments. 

By explicitly going through the six colour configura- 
tions listed above it is easily checked that 



w(x) = — e • Q(x), 
16 



(5.15) 



where the cross-staggered operator p8[ is defined by 

Q(x) = ±Mx) - <r 3 (x)] x [<r 2 (x) - a 4 (x)]. (5.16) 

Since Q(x) is manifestly invariant under 90° rotations of 
the colours around x, Eq. (5.15) is seen to hold true for 
any distribution of the colours around a given vertex. 

In order to find the coarse grained version of io(x) we 
express it as a function of the height field h(x). First 
note that the microscopic operator u>(x) is uniform in 
each of the ideal states of the four colouring model. As 
such it defines a function on the ideal state graph u>(h), 
where h £ T is the coarse grained height. Furthermore, it 
is a periodic function of h and it can therefore be written 
as a Fourier sum: 



w(h) = w e exp(ie • h) 



eeTZ' 



The electric charges appearing in the sum take their val- 
ues in the sub-lattice IZ^ C 1Z* , which is the lattice recip- 
rocal to the lattice of periods of w(h). In the continuum 
limit the coarse-grained height h is promoted into the 
height field h(x), and the scaling limit of the operator w 
is obtained by replacing h by h(x) in Eq. (5.17). There- 
fore u>[h(x)] is a sum of vertex operators, 



Mh(x)] = 



,exp(ie • h(x)) , 



(5.18) 



of which only the most relevant one(s) are kept in the 
effective action. Since the relevance of an operator is 
determined by its scaling dimension we turn to this cal- 
culation next. 



1. Dimensions of charge operators 



In the Coulomb gas formalism operators are associated 
with either electric or magnetic charges. Electric opera- 
tors are vertex operators exp(ie • h) and they appear as 
the scaling limits of microscopic operators in the FPL 2 
model that can be expressed as local functions of the 
colours; the loop-weight operator is one example. 

Magnetic operators on the other hand cannot be ex- 
pressed as local functions of the height but can be 
thought of as a constraint on the height field that gen- 
erates a topological defect of strength m. If x is the 
position of the defect core then the net height increase 
around any loop that encloses x is m (assuming no other 
defects are encircled). Geometrical exponents for loops 
in the FPL 2 model are given by dimensions of electric 
and magnetic operators in the associated Coulomb gas. 

For an operator that has total electromagnetic charge 
(e, m), where e = (ei, e2, 63) and m = (m 1 , to 2 , to 3 ), the 
scaling dimension is the sum of the electric and magnetic 
dimensions, 4 



2x(e, m) = — 



—E a {E a -2E 0a )+g a (M a ) 2 



where 



Ei — ei — e3 , £2 



a 2 



(5.19) 



(5.20) 



(5-17) and 



M 1 = - (to 1 - to 3 ) , M 2 = to 



1 M 3 = -(m 1 +m 3 ) 



(5.21) 



are the electric and magnetic charge vectors in the basis 
in which the elastic term in the action is diagonal. Since 
the magnetic charges are given by height differences they 



must transform according to Eq. (5.7), whilst the electric 



charges transform in a dual fashion (cfr. their appearance 
in the vertex operators). 



4 The derivation of Eq. ( |5.19| ) is an exercise in Gaussian in- 
tegration and is reviewed in Appendix 
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2. Loop ansatz 

With the dimension formula in hand, we can settle 
the issue of the most relevant operators appearing in the 
Fourier expansion of w(h); see Eq. (5.18). There are 



twelve vertex operators to choose from corresponding to 
the twelve (llO)-type vectors in the bcc lattice TZ*; these 
are the shortest vectors in the lattice TZ* W . To find which 
of these electric charges minimise x(e, 0) (Eq. (|5.19| )) it is 
convenient to first consider the simpler case of the FPL 2 
model for n\> — n g . 

For the FPL 2 model with equal fugacities for the black 
and grey loops the effective action is considerably simpli- 
fied. Namely, in this case the cyclic permutation of the 
colours, 



(A,B,C,D)^(B,C,D,A) : 

(zx,z 2 ,z 3 ) <-> (-zi,z 3 ,-z 2 ) 



(5.22) 



does not change the vertex weight A, and is thus an addi- 
tional symmetry of the action S. Thi s sy mmetry implies 
that if 13 — and if 22 = if 11 in Eq. (5.6). Consequently 
there is only one elastic constant, if = if 11. This then 
simplifies the formula for the dimension of an electromag- 
netic charge, 



2x(e, m) 



2nK 



e • (e - 2e ) + —\ 



(5.23) 



where from Eq. ( 5.1 2| ) it follows that the background 
charge in this case has only one non-zero component, 
eo = — 7r(eb, 0, 0). Now it is a simple matter to check that 
of the twelve (llO)-type vectors in the lattice of electric 
charges TZ* , the four charges 

P « = (- 



3 (2) 
3 (3) 

,(4) 



= (-7r,0,+7r), 
(-7r,0,-7r), 

(-7T, +7T,0), 
(-7T, -7T,0) 



(5.24) 



are degenerate in dimension and they minimise 2cc(e,0). 
These are therefore the electric charges of the vertex op- 
erators that are kept in the action. 

Now we turn to the loop ansatz which states that the 
operator u>(h) is exactly marginal in the renormalisation 
group sense. This is the statement that the loop weight 
does not renormalise at large scales. The geometrical 
meaning of this becomes obvious when one realises that 
the number of loops inside a domain of size p, whose 
linear size is comparable to p, is thermodynamically con- 
jugate to the loop weight at scale p. Thus the loop ansatz 
states that the number of large loops does not grow with 
scale (more precisely it is sufficient to assume that it does 
not grow faster than any power of the scale) . The analo- 
gous statement can be proven rigorously for critical per- 
colation where it is the source of hyperscaling p5) . 



The assumption that there is of order one loop at ev- 
ery scale is linked to the variance of the height differ- 
ence between two points in the basal plane, separated by 
a macroscopic distance |x|. Namely, if we assume that 
when going from one point to the other there is of or- 
der one contour loop that is crossed at every scale, and 
further assuming that the directions of these contours 
are independent from scale to scale, it follows from the 
law of large numbers that the variance of the height dif- 
ference grows as the number of contours crossed, that 
is as log(|x|). This of course is nothing but the large 
|x| behaviour of ((H a (x) - H"(0)) 2 ) calculated in the 
Gaussian model of Eq. ( |5.8| ) . 

The loop ansatz, or in other words the marginality hy- 
pothesis for the loop weight operator, simply translates 
into a statement about its scaling dimension: 

ar(eW,0) = 2 



1,2,3,4. 



(5.25) 



This, using the dimension formula Eq. ( 5.23 ), leads to a 
formula for the single elastic constant K. 

In the general case ny, ^ n g , the scaling dimensions of 
the four electric charges identified above are 

1 - e b 



a;(e (3) ,0) 



z(e«. 
x(e^. 
x(e^. 



0) 
0) 
0) 



,9i 



1 



7T 

4 



.93 

1 - 



2e b 



,9i 



1 

.92 



(5.26) 

2e„ 



.9.:! 



the last two remain degenerate in dimension. The di- 
mensions of the first two charges are also equal due to 
the "duality" transformation of the FPL 2 model which 
exchanges the two flavours, n b «-> n g . This transforms 
the microscopic heights zn — > — Z2 and 2:3 — ► —Z3 (and 
similarly for the appropriate components of the height 
field). Furthermore, the elastic constants ifn and if 22 
in Eq. ( |5.6|) ar e unchanged, whilst if 13 — > —if 13. Finally, 
from Eq. ( |5.9| ) it follows that the duality transformation 
exchanges the couplings g\ <-> 33 thus rendering e*- 1 ) and 
e ( 2 ) degenerate in dimension, as the FPL 2 model is self- 
dual. 

Unlike the case of rt b = n g , the loop ansatz in the gen- 
eral case requires that at least two of the electric charges 
eW (i = 1,2,3,4) remain marginal, thus enforcing the 
non-renormalisability of the two fugacities n b and n g . 
If we now further assume that these charges are unre- 
lated by the "duality" transformation described above, 
it follows that in fact all four are marginal. The three 
couplings are then simply calculated by setting the right 
hand sides of Eq. ( 5.26Q equal to 2. We find: 







9i = 


*( 
2 V 






.93 = 


-( 

2 V 


1 


1 


92 


,9i 



eb), 
(l-e g ), 



(5.27) 
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One final comment is in order. The relation l/g-2 = 
l/gi + 1/53 comes as somewhat of a surprise, as it was 
not anticipated on symmetry grounds. Of course, since a 
particular point in the critical region of the FPL 2 model 
is determined by two parameters, and ro g , one rela- 
tion between the three couplings is to be expected. It 
is therefore an interesting open question whether a criti- 
cal loop model can be constructed in which gi would be 
unconstrained. 5 

With the values of the couplings gi,<?2, and 53 in 
hand, as well as the form ula fo r the scaling dimensions of 
charged operators, Eq. ( 5.19 ), we are fully equipped to 
calculate critical exponents of the FPL 2 model. In par- 
ticular, in the next section we calculate the formulae for 
the central charge and the geometrical exponents associ- 
ated with loops as a function of the loop fugacities, n\, 
and rig, for the whole critical region of the model. 



c = 3-6 



1 — e b 1 — e„ 



(6.4) 



where we recall that — 2 cosset,) and similarly for n g . 
In Table | the numerically calculated values of the con- 
formal charge are compared to the above formula, and 
excellent agreement is found. 



VII. GEOMETRICAL SCALING DIMENSIONS 



VI. CENTRAL CHARGE 

We now turn to the calculation of the central charge 
in the critical region, < n-b,n g < 2. Exactly at the 
point (rib,n g ) = (2,2) the background charge vanishes, 
eo = 0, and the action consists only of the elastic term 
<Se given by Eq. ( |5.8| ) . Since this is then simply a theory 
of three free massless bosonic fields we conclude that, in 
this case, c = 3 p8[ . 

For a general value of the background charge this gen- 
eralises to 



c = 3 + 12x(e ,0) 



(6.1) 



One way to rationalise the factor of 12 is to compare the 
coefficients of the finite-size corrections in the well-known 
formulae |37|,i38| 



/o(oo) 



ML) - fo(L) = 



(6.2) 
(6.3) 



where fo,i(L) is the free energy density on a cylinder of 
circumference L, the subscript referring to the vacuum 
and i to the case when an operator of scalin g dim ension 
Xi is inserted. The physical meaning of Eq. ( |6.l[ ) is that 
the presence of the background charge — heo and — eo at 
the two ends of the cylinder - lowers the free energy and 
with it the central charge. 



Now using the dimension formula, Eq. (5.19), and in- 
serting the values of the couplings g a from Eq. (5.27), we 
arrive at 



A. Two-string dimension 



In addition to the central charge, the Coulomb gas rep- 
resentation of th e lo op model provided by the Liouville 
field theory, Eq. (5.2), allows for the evaluation of various 
geometrical scaling dimensions. As an example of such a 
quantity, consider the probability G2(r) that two points 
separated by a distance r lie on the same, say, black loop. 
In the critical phase we expect this probability to decay 
as G2{r) ~ r ~ 2x2 , which defines the scaling dimension 
X2 ■ Since a black loop is represented as a sequence of 
alternating A and B-coloured edges it follows from the 
colouring constraint that the microscopic heights z just 
outside this loop differ by integer multiples of C and D 
only. In other words, a black loop is a contour loop for 
the component of the height along the direction perpen- 
dicular to both C and D, i.e., the (1,0,-1) direction in 
height space. Similarly the grey loops are contour loops 
for the height component along the (1, 0, 1) direction. 

It has been argued that the scaling dimension govern- 
ing the probability that two points belong to the same 
contour loop of a random Gaussian surface equals 1/2, in- 
dependent of the stiffness |3{|. Thus, for (nb,n g ) = (2, 2) 
when eo = and the effective field theory is Gaussian, we 
expect X2 — 1/2. For other values of the fugacities the 
Gaussian theory is modified by the background charge 
and the same argument cannot be made. 



5 This possibility was suggested to us by D. Huse. 
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a) 




b) 




FIG. 7. Defect configurations used to calculate geometrical 
exponents xi (a) and X2 (b) in the FPL 2 model. In (a) there 
is a single oriented black loop segment and a single oriented 
grey loop segment propagating from to r, whilst in (b) there 
are two oriented black loop segments between and r. 

A more illuminating way of making contact with the 
interface representation is to view G2(r) as a two-string 
correlation function associated with defect configurations 
where two black strings emanating from the origin an- 
nihilate one another at a distant point r; see Fig. f?]b. 
This can be accomplished by rewriting G2{r) as Z(r)/Z, 
where Z is the partition function defined by Eq. ( |4.3| ), 
and Z(r) is similarly defined but with the summation 
restricted to those configurations Q' r where an oriented 
black loop passes through the points and r. Now con- 
sider reversing the direction of one half of the loop, so 
that instead of having one oriented loop passing through 
and r we have two oriented loop segments directed from 
to r |J. This corresponds to the introduction of de- 
fect configurations at these two points, where we have 
violated the edge-colouring constraint. At we find a 
(C, D, A, A) configuration of colours which in the height 
language corresponds to a vortex of strength 



m 2 



A-B 



-2,0,2) . 



(7.1) 



The strength of the vortex (its Burgers charge) at is 
calculated as the total height change around 0. Similarly, 
at r we have the corresponding antivortex (B,B,C,D) 
of strength m 2 as illustrated in Fig. ^|b. 

In order to calculate Xi for general values of the loop 
fugacities we have to take into account the effect of the 
complex phase factors associated with oriented loops. 
Namely, when one or more, say, black strings are as- 
sociated with a vortex-antivortex configuration, spurious 
phase factors exp(±«7reb) will arise whenever a black loop 
segment winds around one of the vortex cores [pi ; for ex- 
ample, in Fig. [?]b one of the two black strings winds once 
arround point r. The spurious winding phase can be re- 
moved by inserting the vertex operator exp(ieb • h) at 
the positions of both vortex cores. Since a black loop has 



alternating A and B colours the electric charge eb must 
satisfy 



e b • A = 7re b , 
e b • C = 0, 



e b • B = -7re b , 
e b • D = . 



(7.2) 



Similarly, if there are grey strings propagating between 
two vertices the spurious phase factors associated with 
winding configurations are corrected with vertex opera- 
tors whose electric charge e g is determined by 



0. 



D 



B 



-0, 



(7.3) 



Using Eq. (4.2) for the colour-vectors we find, 



eb 



-(e b ,0, -e b ), 



: (e g ,0,e g ). (7.4) 



Going back to the two-string operator we conclude that 
it has total electromagnetic charge (eb,m2). 

Finally, from the general expression for the scaling di- 



mension of an electro-magnetic operator, Eq. (5. IE), it 
follows that 



2x 2 = 2x(e h , m 2 ) = (1 - e b ) - 



1 - e b 



(7.5) 



In Table |V| exact values of X2 calculated from this for- 
mula are compared to numerical results, and excellent 
agreement is found. 

Interestingly the expression for X2 is independent of e g , 
i.e., it is not affected by the fugacity of grey loops. This 
observation conforms to our understanding of the scaling 
of compact polymers. The compact polymer problem is 
recovered in the limit nb — > in which case there is a 
single black loop on the lattice. Since the loop fills space 
its Hausdorff dimension is necessarily D — 2. Scaling 
tells us that EfJ 



D 



X2 



(7.6) 



from which the result X2 = follows, independent of the 
fugacity of grey loops. The fact that our formula repro- 
duces this simple result in the = (eb = 1/2) case 
provides a non-trivial check on its validity. 



B. One-string dimension 

The scaling dimension x\ corresponding to one black 
and one grey string propagating between two points on 
the lattice, can be computed in a way that is com- 
pletely analogous to the case of two black strings dis- 
cussed above. (Note that the fully packing constraint 
ensures that if there is a single black string between two 
points then these points are also connected by a grey 
string; see Fig. [7|a.) Choosing one point on the even sub- 
lattice and the other on the odd, leads to the appearance 
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of the defect configuration (A, C, C, D) on both sites of 
the square lattice. These in turn correspond to vortices in 
the height representation with topological charges ±mi, 
where 



m! = C-B = (-2,-2,0) 



(7.7) 



Since strings of both flavours are now present the com- 
pensating electric charge is eb + e g = eo. Hence 



2xi = 2x(eo, mi) 



3 [(l-e b ) + (l-e g )] 

(l-e b )(l-e 8 ) 
(1 - e b ) + (1 - e g ) 



(7.8) 



1 - e b 



1 



There are of course several different ways of choosing the 
defect configurations (in this case, eight), but it should 
hardly come as a surprise that they all lead to the same 
expression for the scaling dimension. 

Unlike X2, x\ depends on both loop fugacities. Going 
back to our original motivation, the compact polymer 
problem (n b = =>■ e b = 1/2), x\ determines the value 
of the conformational exponent 7 = 1 — x%, which de- 
scribes the scaling of the number of compact polymers 
with the number of monomers. We see that depending 
on e g there will be a continuum of 7's. How do we inter- 
pret this? 

First note that the problem of counting the number 
of conformations of a single compact polymer is the case 
rig = 1 (e g = 1/3) which simply assigns equal weights 
to all conformations. Using Eq. ( |7.8| ) this choice leads to 
x% = —5/112 and to the result 7 = 117/112 advertised 
in the abstract. Changing n g (e g ) away from n g = 1, 
on the other hand, has the effect of favouring certain 
compact polymer conformations over others depending 
on the number of loops formed by the uncovered (grey) 
bonds. In this sense the weight assigned to grey loops can 
be thought of as an interaction between the monomers 
of the compact polymer, albeit a peculiar non-local one. 
Nonetheless, it is interesting that this interaction changes 
the scaling properties of the compact polymer leading to 
a continuously varying exponent 7 (more on this in the 
Discussion). 



C. Many-string dimensions 



the defect configurations obtained by violations of the 
edge colouring constraint must necessarily give rise to 
an even number of strings we will only consider the case 
when s b + s g is even. 

Consider first the case s b = 2fc b and s g = 2k g . The 
appropriate magnetic charge is obtained by combining fc b 
vortices with charge A — B = (—2, 0, 2), and k g vortices 
with charge C — D = (—2, 0, —2). The defect with charge 
A— B acts as a source of two black segments, whilst C— D 
is associated with two grey loop segments. We also need 
to introduce the electric charge e b + e g to compensate for 
the extra winding phase associated with the black and 
grey loop segments. The total electromagnetic charge is 
therefore 

[e2fe b ,2fc s , Hi 2 fc b ,2fcJ = (7.9) 
[eb(l - <5fc b ,o) + e g (l - 4 g , ),-2(fc b + fe g ,0, k g - k b )], 

and from the dimension formula, Eq. ( [5.19] ), we find 



(1 - e h )kl + (1 - e g )fc 2 + 



1 - eb 



(1 



°k h ,0 



1 



e„- 



-(i-4 g)0 ). 



(7.10) 



This formula generalises Eq. (7.5). 

Similarly, for s b = 2fc b — 1 and s g = 2k g 
tromagnetic charge is 

[e2fe b -i,2fe g -i 5 ni2fe b _i ) 2fe g -i] = 
[e , -2(fcb + fcg - 1, 1, k g ~ k h )}; 



1 the elec- 



(7.11) 



the magnetic charge is obtained by combining fc b — 1 de- 
fects of charge A — B, k g — 1 defects of charge C — D, 
and a single defect of charge C — B which produces the 
remaining single black and grey strings originating from 
the same vertex. The scaling dimension is found to be 

2x2fc b -l,2fc g -l = 

J [(1 - e b )(2fc b - l) 2 + (1 - e g )(2fc g - l) 2 ] + (7.12) 
(l-e b )(l-e g ) 



(1 - e b ) + (1 - e g ) 



1 



eb 



1 — e„ 



This generalises the e xpre ssion given in Rcf. 
rectly reduces to Eq. (7.8) for & b ,/c g = 1. 



and cor- 



Thc dimensions x\ and X2 given above are contained 
in a more general set of string dimensions x SbjSg govern- 
ing the probability of having s b black loop segments and 
s g grey loop segments propagating between two points 
on the lattice [H. More precisely, we consider two mi- 
croscopic regions centred around points separated by a 
macroscopic distance, one region being the source and 
the other the sink of the oriented loop segments. Since 



D. Thermal dimension 

We now turn our attention to the thermal scaling di- 
mension. The FPL 2 model can be thought of as the 
zero-temperature limit of a more general model where we 
allow for thermal excitations that violate the close pack- 
ing constraint. In this sense the temperature variable 
is thermodynamically conjugate to the constraint that 
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every vertex be visited by (say) a black loop. An appro- 
priate defect configuration for computing xt within the 
FPL 2 model is therefore (C,D,C,D). This is a vortex 
of strength 



mr = 2(C + D) = (0,-4,0), 



(7.13) 



and since no strings terminating in the bulk are generated 
there is no compensating electric charge. The scaling di- 
mension is then 



2x T = 2x(0,m T ) = 4 



(l-e b )(l-e g ) 
(1 - eb) + (1 - e g ) 



(7.14) 



The exact values of xt quoted in Table [fj are calculated 
using this formula. The numerical results are in excellent 
agreement. 



E. Boundary-string dimensions 

The simplest example of a string operator that cannot 
be accessed within the formalism presented above is that 
of one black and no grey strings propagating between two 
vertices of the square lattice. Since this configuration has 
an odd number of strings connecting two sites of the lat- 
tice these two sites necessarily reside on the boundary. 

If we define the FPL 2 model on the cylinder, as will be 
the case when we construct its transfer matrix in Sec. IX , 
a single black string can be enforced to run along the 
length of the cylinder if its circumference is chosen odd. 
Taking our cue from the formulae derived above for the 
bulk string operators we guess the formula 



X = - 

8 



1 , l-e b 



1 



2 1 - e b 



(7.15) 



III. X is the 



from the numerical results shown in Table 
scaling dimension of the boundary operator which corre- 
sponds to a single black (or grey) string. 

The Coulomb ga s int erpretation of the second and 
third term in Eq. ( 7.15) is rather apparent when one 



compares them to Eq. (7.5). The second term can be 
rationalised as coming from a mag netic charge (—1, 0, 1) 
which is half the charge ni2 in Eq. (7.1), associated with 
two black strings; this is saying that we have a partial 
dislocation generated at the boundary. The third term 
is due to the compensating electric charge eb for a single 
black string, same as in the two-string case. 

The first, constant term does not have an immediate 
interpretation. A possible scenario is that it is due to the 
boundary condition imposed on the height by virtue of 
having a cylinder of odd circumference. Namely, a trans- 
lation along the periodic coordinate by an amount equal 
to the circumference (L) exchanges an even site for an 
odd site (and vice versa) resulting in a transformation of 
the height: h(x 1 ,a; 2 ) = Ph(x 1 +L,x 2 ). Since P 2 = 1 this 
boundary condition can be thought of as an insertion of 



a twist operator into the partition function. The twist 
operator has dimension 1/8 regardless of the stiffness of 
the interface H |. 

The above considerations permit us to calculate the 
scaling dimension for the general case of an odd num- 
ber of strings. For definiteness we consider the case of 
,s b = 2fc b ~ 1 and s g = 2fc g . The magnetic charge per- 
taining to this situation is found by combining 2/c b — 1 
defects of charge ^(A — B) with k g defects of charge 
C D, totaling 

[ e 2fc b -l,2fc g , ni2fc b -l,2fcj = (7-16) 

[e b + e g (l - 4 g ,o), (1 - 2fc b - 2fe g , 0, 2fc b - 2fc g - 1)]. 

Taking into account the contribution from the twist 
operator, i.e., adding 1/8 to the result obtained from 
Eq. ( 5.19| ), the scaling dimension is then 



2x 



2fc b -l,2fc„ 



i + i[(l-e b )(2fc b -l) 2 



(l-e g )(2fc g ) 5 



(7.17) 



■0--Sk.,o) 



l-e b 1 



F. Complete spectrum of string dimensions 



Finally, the results of Eqs. ( ]7T0| ), fTVZ} ) and ( fTU^ ) 
can be combined into a single equation for the scaling di- 
mension of a string operator that corresponds to Sb black 
loop segments and s g grey loop segments: 



2.E 



Sb,S 

1 



+ 4 [(1 - *0-2 + (1 - e g )4] 



(7.18) 



l-eb 



(l-*- b ,o) + 



-(l-f-.,o) 



A (2) ,(2) (1 - e b )(l - e g ) 



(l-e b ) + (l-e g )' 



(2) 

here we have defined o\ ? = o i=J( - mod 2 )- 



VIII. TERMINATION OF CRITICAL 
BEHAVIOUR 

In the preceding sections we have developed an effec- 
tive description of the critical phase of the FPL 2 model in 
the form of a field theory. This theory has to break down 
at large values of the loop fugacity since in this case a 
typical state of the model will consist of small loops only, 
i.e., a power-law distribution of loop sizes will be absent. 
That this indeed happens can be seen from the Liouville 
field theory itself as it carries the seeds of its own demise. 
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The mapping of the loop model to an oriented loop 
model for , n g < 2 works equally well for > 2 or 
n g > 2. From Eq. (4.1) it follows that in the latter case 



at least one of the parameters, eb or e g , will be pure 
imaginary. This affects the Liouville potential which for 
«b > 2 or n g > 2 becomes a relevant perturbation to the 
(modified) Gaussian action Se + £b . 

To understand how this comes about we consider the 
simple case provided by the = n g FPL 2 model, dis- 



cussed in Sec. V C 2. Namely, as we increase the value of 



the loop fugacity we expect small loops to be favoured 
and the stiffness K of the interface to grow. In the 
critical phase this is offset by the decrease in the back- 
ground charge in a way that leaves the Liouville poten- 
tial marginal. Now when the loop fugacity exceeds 2 the 
background charge eo = — 7r(eb, 0, 0) becomes pure imag- 
inary and the dimension of the Liouville potential 



XL= 2 



7r 1 - e b 



K 



.1 



can no longer stay marginal; here xl = x(e^\0), where 
the charges e^' are given in Eq. (5.24), and their scal- 



ing dimensions are calculated from Eq. ( 5.23| ). In fact, 
assuming that the stiffness K continues to increase with 
the loop fugacity for rib = n g > 2, xl turns complex 
with a real part that is smaller than two, rendering the 
Liouville potential relevant. 

If we make the usual assumption of no intervening fixed 
points, the relevant Liouville potential will generate a fi- 
nite correlation length and the loop model will no longer 
be critical. The correlation length has the physical in- 
terpretation of the average size of a loop in the system. 
This scenario has been confirmed for the fully packed loop 
model on the honeycomb lattice, from the Bethe ansatz 
solution of this model jl2| . 

A different view of the non-critical region of the FPL 2 
model is provided by the locking potential V(h). Namely, 
the discrete nature of the microscopic heights can be 
taken into account in the field theory by a negative po- 
tential in height space that is peaked around the flat, 
ideal states. As such, this potential is uniform on the 
ideal state graph and can therefore be expanded in a 
Fourier series. Examination of the most relevant ver- 
tex operators in this series p8[ reveals that they are the 
same as the ones for the loop- weight (Liouville) potential, 
w(h). Therefore, just like w(h), the locking potential in 
the non-critical region of the phase diagram is a relevant 
perturbation. Thus, it will lock the interface in one of 
the ideal, flat states. In this flat phase the height fluctu- 
ations are bounded (as opposed to being logarithmically 
divergent) which is just another way of saying that large 
contour loops are exponentially suppressed. On the other 
hand, in the critical region of the FPL 2 model the lock- 
ing potential is marginal as it would be for an interface 
model at the roughening transition Kj| . This might in- 
dicate that the whole critical region of the FPL 2 model 



can be understood as a manifold of essential singulari- 
ties in some more general model, as was the case for the 
honeycomb FPL model QH . 



As advertised in Sees. VI and VII our results for the 



central charge and a number of the geometrical scaling 
dimensions have been very accurately confirmed by trans- 
fer matrix calculations. Before turning to a discussion of 
our numerical results we describe the particular repre- 
sentation of the transfer matrix used to obtain them. 



IX. CONSTRUCTION OF THE TRANSFER 
MATRIX 

To construct the transfer matrix for the FPL 2 model 
on a cylinder of circumference L we write the partition 
function as 



Z (M) = 



E 



(9.1) 



where the length of the cylinder M has been explicitly 
indicated. Periodic boundary conditions are imposed in 
the horizontal direction, whereas the bottom and the top 
row of the cylinder have open boundary conditions and 
hence terminate in L dangling edges. We recall that the 
restriction of the summation to the set of fully packed 
graph configurations Qm implies that locally the vertices 
are constrained to have one of the six appearances shown 
in Fig. ||. In the first four possible vertices the loop seg- 
ments do not cross, whilst in the last two vertices the two 
flavours intersect. The global constraint that all loops be 
closed in the limit of an infinite system means that loop 
segments cannot terminate in the bulk but only at the 
dangling edges in the top and bottom rows. 

A typical loop configuration for a cylinder with L = 6 
and M — 12 is shown in Fig. ||. The horizontal num- 
bering pertains to the vertices, whilst in the vertical di- 
rection it is more convenient to label each row by the 
number of the vertex immediately below it. Accordingly 
the labels and M refer to the bottom and the top row 
of dangling edges respectively. We shall soon see that the 
inclusion in Qm of one or more strings running between 
the dangling edges of row and M helps us access the 
geometrical exponents of the model. In particular, the 
configuration of Fig. || having one such string of each 
flavour furnishes a contribution to the scaling dimen- 
sion xi which determines the conformational exponent 

J=l-Xi. 



A. Connectivity basis 

The construction of a transfer matrix (TM) for 
Eq. ( |9.lD appears to be obstructed by the non-locality 
of Ni (i = b, g). The key to solving this problem is to 



1G 



write the TMs in a basis of connectivity states compris- 
ing information about how the dangling ends of row M 
are pairwise interconnected in the preceding rows and, 
if strings are present, information about the positions of 
such strings. In addition the connectivity states must 
keep track of the particular flavour of any loop or string 
segment terminating in row M. Our construction gener- 
alises the work of Bldtc and Nightingale for the Q-state 
Potts model |45|,E6] and that of Blote and Nicnhuis for 
the O(n) model 47 to take the extra flavour information 
into account, and our notation is consistent with that of 
these authors. 
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It is essential to be able to represent a given connec- 
tivity state both in an index representation giving direct 
access to the flavour and connectedness information just 
mentioned, and in a number representation assigning an 
integer in the range 1,2,..., C , ^ St " Sg - ) to the state under 
consideration. The latter representation enables us to 
enumerate the entries of the TM, whilst the former al- 
lows us to determine the number of loop closures when 
going from one connectivity state to another and hence 
the value of a particular entry in the TM. Here C^ t " Ss " > is 
the number of distinct connectivity states for a cylinder 
of width L accommodating S{ strings of flavour i = b, g. 
The construction of these two representations, the map- 
ping between them, and the evaluation of the C'£ ib ' Sg - ) 
for (s b ,s g ) = (0,0), (1,0), (1,1) and (2,0) is deferred to 
Appendix |b|. 



FIG. 8. A typical loop configuration for L = 6 and M = 12. 
The dashed lines along the left and the right boundaries il- 
lustrate the periodic boundary conditions. Horizontally the 
vertices are numbered from 1 to L, whilst vertically the rows 
are labeled by the number of the vertex immediately below 
them. This particular configuration is constrained to having 
precisely one string of each flavour spanning the length of the 
cylinder, and hence it contributes to the geometrical expo- 
nent x\. To the right we show the index representation of 
the connectivity state pertaining to each row (see Appendix 
[b] for details). Any valid configuration can be interpreted as 
a "jigsaw puzzle" assembled from the six "pieces" shown in 
Fig. H. Note that when laying down the first row of this puz- 
zle it must be stipulated how the dangling edges of row 0, 
which are not part of a string, are pairwise interconnected 
below that row. These implicit connections as well as their 
counterparts in row M have been depicted by dashed loop 
segments. 

Designating the connectivity states by Greek letters 
we can write the partition function as a sum of restricted 
partition functions 

Z(M) = £ Z (M) = £ £ ^ ^ { g M))n N bn N^ (9 2) 

P 0m 

where 4>{Gm) is the connectivity of the L dangling edges 
of row M, and 5(i,j) is the Kronecker delta. Now con- 
sider appending another row to the cylinder, giving us 
a total graph configuration Gm+i = Gm U G' ■ Evidently 
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the connectivity of the dangling edges of row M + 1 is 
determined solely by that of the preceding row and by 
the appended subgraph Q' 



<KGm+i) = i/>(<KGm),g'). 



(9.3) 



Letting N- denote the number of loop closures induced 
by Q' we arrive at the relation 



Z (M + i) = £ 6 (a,<f>(G M+ i)) 

Gm + i 

= EE*(m&*)) 

5(a^{ct>{GM),G')) 

Q'\Gm 

Err 7 {M) 



where the transfer matrix is defined by 



N b +N' N e +N' 



hi-". 



K K 



Q'\Qm 



iw K K 



(9.4) 



(9.5) 



The notation Q'\Qm means that the summation is con- 
strained to those subgraphs Q' that fit the dangling edges 
of Gm- 



B. Single- vertex decomposition 



A quintessential step in the practical implementation 
of the TM is its decomposition into matrices each corre- 
sponding to the addition of a single vertex, 



-L-1 



(9.6) 



Here the single- vertex matrix Tj, which adds the vertex 
at horizontal position i of the new row, has the advan- 
tage of being sparse, and we shall soon see that it has at 
most three non-zero entries per column. This property 
leads to a dramatical reduction of the time and storage 
requirements for the calculations. 

As was the case in the O(n) model E?]], a minor com- 
plication arises due to the fact that the addition of the 
first vertex of a new row increases the number of dangling 
edges from L to L + 2. This is illustrated in the left part 
of Fig. |^. Upon addition of further vertices the number 
of dangling edges is kept fixed at L + 2, until the L'th 
vertex completes the row, and we are back at L dangling 
edges. Thus the dimensions of the single-vertex matrices 
are C L+2 x C L for T 1; C L+2 x C L+2 for T 2 , . . . , Tl-i, 
and C L x C L+2 for T L . 
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FIG. 9. Adding the first vertex of the (M + l)'th row increases the number of dangling edges from L to L + 2. The labeling 
of the "active" edges (filled circles) before and after addition of the new vertex (shaded) is as shown in the left part of the 
figure. The part of the lattice relevant for determining which of the vertices of Fig. [] fit onto a given connectivity of row M, 
has been depicted in the right part of the figure. This information constitutes the vertex rules, and is explained in the text. 



In Fig. [9] we illustrate the action of Ti on z( M ' in de- 
tail. To ensure that row M + 1, when completed, will 
have the same labels on its dangling edges as was the 
case in the preceding row, the solid dots illustrating the 
"active" dangling edges must be relabeled as shown in 
the lower left part of the figure. Shown to the right are 
the three possible choices of vertices fitting onto a black 



loop segment terminating at the dangling end 1 of Z^ M > . 
There are thus three non-zero entries in each column of 
Ti. Since no loop closures of either flavour can be in- 
duced (N^ — iVg = in Eq. (9.5)) all these entries are 
unity. Similar considerations hold true when the loop 
segment to be fitted is grey, and the vertex rules can be 
read off from the figure by interchanging the two flavours. 



18 



□ 



FIG. 11. Completing the (M + l)'th row by adding the L'th vertex (shaded) brings the number of dangling edges back to 
L. The labeling is now consistent with that of the preceding rows. Vertex rules, shown to the right, now include a disallowed 
configuration. Namely, in situation (a) none of the vertices of Fig. |^ fit in, and the corresponding entry of the transfer matrix 
must be forced to zero. Situations (b), (c) and (d) offer various possibilities for a black loop closure. 
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When strings are present a few modifications of the 
above rules are necessary. In situation (a), if one of i\ 
and ig equals -1 and the other is positive, the left-over 
partner to the non-string black segment must be made 
the new string. And if both i| and i\ equal -1 the cor- 
responding entry of T2 must be forced to zero, since two 
strings cannot be allowed to annihilate. 

Finally, consider closing the (M + l)'th row through 
the action of T^, as depicted in Fig. [ll]. The labels L+ 1 
and L + 2 now disappear, and as far as the labeling goes 
the system is back in its original state. Each column of 
Tl has at most one non-zero entry per column, as wit- 
nessed by the vertex rules displayed in the right part of 
the figure. Once again, only half of the vertex rules are 
shown, and the other half is found by interchanging the 
two flavours. 

In situation (a) no vertex of Fig. || can fit onto the 
three dangling edges at positions L, L + 1 and L + 2. 
The corresponding entry of must therefore be forced 
to zero. Situations (b), (c) and (d) leave us to determine 
whether, for a given connectivity of the L + 2 dangling 
edges, a black loop closure occurs or not. The handling in 
terms of the index representation is exactly as described 
above. 



Z (M) 



^ T (0,0) 



M 



(10.2) 



It is well-known that conformal invariance relates the 
amplitude of the 1/L 2 corrections to /o°'°' ) (oo) to the cen- 
tral charge c |37|] . A further (non-universal) 1 /i 4 correc- 
tion due to the operator TT, where T denotes the stress 
tensor, must also be present in any conformally invariant 
syste m [H . It is therefore found in a number of cases 
||,|9|j50j that fits of the form 



ft°\L) = fr> (00) 



(0,0), 



7TC A 

6I2 + U 



(10.3) 



yield very rapidly converging estimates for c. An effi- 
cient application of Eq. (10.3) is to determine c from 
parabolic least-squares fits of the finite-size data against 
1/L 2 Ha 

In Table Q the results of such fits including the data 
points for L < L < L max are shown as a function of 
Lq. Numerically we were able to access L max = 14, in 
which case the largest single- vertex TMs have dimension 
~ 7 • 10 6 (see Table VII). The extrapolation of the esti- 
mants c(Lq, L mayi ) to the limit of infinite Lq is assumed 
to take the form of a power law 



X. NUMERICAL RESULTS 

A. Central charge 

The reduced free energy per vertex in the limit M — » 00 
of an infinitely long cylinder is given by 



/ (0 " 0) (i) = , lim 



1 



M^oo LM 



lnTr^W 



1 
L 



InA 



(0,0) 
' 



(10.1) 



where A Sb,Ss ^ is the largest eigenvalue of T^ b,Sg ^. The 
partition function for a cylinder of length M is found by 
iterating the no-string TM 



c(L ,L max ) = c + kL p , 



(10.4) 



at least within an asymptotic regime of large enough Lq. 
As is evident from Table | the last three estimants usu- 
ally exhibit monotonicity, thus allowing us to fix the con- 
stants c, k and p. When this was not the case, or when- 
ever the power p thus obtained was to small to produce a 
reliable extrapolation the Ising-like value p = 2 was used 
by default to extrapolate the last two estimants. An error 
bar for this type of fit can be estimated from the varia- 
tion among the individual estimants. The extrapolants 
are invariably in excellent agreement with our analytical 
results, the relative deviation being typically of the order 
10~ 3 . 



rib 


n g 


c(4, 14) 


c(6, 14) 


c(8, 14) 


c(10, 14) 


Extrapolation 


Exact 


0.0 


0.0 


-2.8943 


-2.8861 


-2.9220 


-2.9514 


-3.0037 


-3.0000 


0.5 


0.0 


-1.8528 


-1.7641 


-1.7716 


-1.7873 


-1.8152 


-1.8197 


0.5 


0.5 


-0.7295 


-0.6249 


-0.6159 


-0.6220 


-0.6328 


-0.6395 


1.0 


0.0 


-1.0012 


-0.9542 


-0.9636 


-0.9761 


-0.9983 


-1.0000 


1.0 


0.5 


0.1341 


0.1877 


0.1924 


0.1895 


0.1843 


0.1803 


1.0 


1.0 


0.9918 


0.9969 


0.9986 


0.9999 


1.0004 


1.0000 


1.5 


0.0 


-0.3765 


-0.3669 


-0.3817 


-0.3923 


-0.4111 


-0.4124 


1.5 


0.5 


0.7652 


0.7746 


0.7729 


0.7715 


0.7690 


0.7678 


1.5 


1.0 


1.6215 


1.5818 


1.5778 


1.5806 


1.5856 


1.5876 


1.5 


1.5 


2.2541 


2.1691 


2.1581 


2.1627 


2.1709 


2.1751 


2.0 


0.0 


0.0706 


0.0549 


0.0342 


0.0235 


-0.0019 


0.0000 


2.0 


0.5 


1.2209 


1.1937 


1.1868 


1.1861 


1.1849 


1.1803 


2.0 


1.0 


2.0792 


2.0002 


1.9899 


1.9937 


2.0005 


2.0000 


2.0 


1.5 


2.7139 


2.5919 


2.5737 


2.5781 


2.5859 


2.5876 
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2.0 2.0 3.1629 3.0121 2.9885 2.9936 3.0027 3.0000 

TABLE I. Estimants c(Lo, L max ) for the central charge are obtained from parabolic least-squares fits against 1/L 2 using the 
numerical data for Lq < L < L max . The extrapolation in Lq is described in the text. 



The results for c are shown for all integer and half- 
integer values of S [0,2]. Because of t he symmetric 
appearance of the two flavours in Eq. (9.1) only rib > n g 
need be considered. For either rib = 1 or n g = 1 the 
FPL 2 model reduces to the simpler FPL model earlier 
considered by Batchelor et al. fbifl , and for rib = n g we 
recover another special case recently investigated by one 
of us fl. 



B. Thermal scaling dimension 

A further prediction of conformal invariance is that the 
finite-size scaling of the first gap in the eigenvalue spec- 
trum of T^ ' ) is related to the thermal scaling dimension 



f[ o ' a) (L)-fr>(L) 



(0,0) 



where f[ ' ^ is found from the next-largest eigenvalue of 
T(°<°) through /f' 05 = -±lnAf' 0) . These computations 



2ttxt 



(10.5) 



were also carried through for even L up to L max = 14. 
In this case as well the convergence of the estimants can 
be c onsid erably sped up by including a 1/L 4 term in 
Eq. ( 10.5 ) and performing parabolic least-squares fits ver- 
sus 1/L 2 . 

The results for xt as displayed in Table [0] again agree 
with those of the previously studied special cases [Qff) ■ 
The data for (nb,n g ) = (0,0) merit a special comment. 



Monitoring the three leading eigenvalues A 
Ao 0,0 ' as a function of n for — r 



(0,0) A (o,o) 



and 



n we found that 



A^ 0,0 "* and X^'"' are exactly degenerate for all n down to 



v (o,o) 



n ~ 0.20. Hereafter A l j '°' ) splits off from X^'"' and even- 
tually becomes degenerate with A °'°' ) at n = 0. Because 
of this level crossing it thus seems very likely that near 
(m, 712) = (0, 0) the thermal eigenvalue should be related 
to the gap /2°'°' ) (L) — /o°'°' ) (L) . Comparison with the ex- 
actly known result xt = 1/2 confirms this suspicion. 
A similar comment holds true near (rib, n s ) = (2, 2), and 
again we find fair agreement with the exact result if we 



(0,0) 



apply Eq. (10.5) to X^'°\ and not to Ap u; (which in this 



(0,0) 



case becomes two- fold degenerate). 



rib 


n s 


x T (4, 14) 


kt(6, 14) 


x t (8,14) 


x T (W, 14) 


Extrapolation 


Ref. [14] 


Exact 


0.0 


0.0 


0.5712 


0.5280 


0.5121 


0.5060 


0.4987 




0.5000 


0.5 


0.0 


0.5704 


0.5535 


0.5452 


0.5417 


0.5366 




0.5372 


0.5 


0.5 


0.5916 


0.5882 


0.5845 


0.5825 


0.5789 




0.5804 


1.0 


0.0 


0.5826 


0.5798 


0.5765 


0.5748 


0.5708 


0.573 (1) 


0.5714 


1.0 


0.5 


0.6204 


0.6227 


0.6218 


0.6211 


0.6199 


0.6200 (5) 


0.6206 


1.0 


1.0 


0.66368 


0.66600 


0.66642 


0.66654 


0.66663 


0.6666 (1) 


0.66667 


1.5 


0.0 


0.5965 


0.6053 


0.6060 


0.6058 


0.6054 




0.6063 


1.5 


0.5 


0.6493 


0.6559 


0.6574 


0.6578 


0.6585 




0.6619 


1.5 


1.0 


0.7782 


0.7094 


0.7108 


0.7115 


0.7130 


0.713 (1) 


0.7146 


1.5 


1.5 


0.8950 


0.7657 


0.7674 


0.7684 


0.7702 




0.7699 


2.0 


0.0 


0.6167 


0.6295 


0.6338 


0.6349 


0.6356 




0.6667 


2.0 


0.5 


0.7481 


0.6878 


0.6913 


0.6927 


0.6945 




0.7345 


2.0 


1.0 


0.8741 


0.7566 


0.7552 


0.7565 


0.7588 


0.76 (1) 


0.8000 


2.0 


1.5 


0.9436 


0.8755 


0.8284 


0.8303 


0.8337 




0.8702 


2.0 


2.0 


0.9996 


0.9850 


0.9400 


0.9200 


0.8876 




1.0000 


TABLE II. 


The thermal scaling 


; dimension xt- 


The extrapolation of the estimants xt(Lq, L max ' 


is described in 


the text. For 



comparison we also show the numerical data for the case of either rib or n g being unity [14]. Due to level crossing the values 
of xt for (rib, n g ) = (0,0) and (2,2) are found from the gap /2°'°'(I/) — fo°'°\L) rather than from f[°'°\L) — /q°'°'(L). 



For rib < 2 the extrapolants are again in excellent 
(~ 10~ 3 or better) agreement with our analytical re- 
sults. For rib = 2 the slower convergence can be at- 
tributed to logarithmic corrections plj ] arising from an 
enhanced number of marginal vertex operators. Indeed, 



of the twelve vertex operators corresponding to the short- 
est vectors in 7£* , Eq. ( 5.18 ), seven stay marginal when 
either rib < 2 or n g < 2. In the general case, when both 
rib < 2 and n g < 2, there are only four marginal vertex 
operators; this is the loop ansatz, Eq. (5.25). 
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C. Dimensions of string operators 

We now turn our attention to the determination of the 
scaling dimensions associated with one or more strings 
spanning the length of the cylinder. The presence of one 
black string corresponds to a height mismatch in the ideal 
states, and the relevant scaling dimension X is therefore 
that of a twist-like operator [|4l|. We have calculated the 
leading eigenvalue of T^ 1,0 ) for odd system sizes up to 
L max = 13 and determined the corresponding estimants 
l(Lg, L max ) by the usual parabolic fits to f^ ,0 \L), 



Coddl 

cfr 
by 



Eq. ( |10.3| ). Estimants X(L Q , L max ) are then defined 



C o dd(i0j imax) 

12 



(10.6) 



wher e the factor of 12 o riginates from a comparison of 
Eq. ( |10.3] ) with Eq. ( 10.5 ). For the central charge c of a n 
even-sized system we use our analytical results, Eq. ( |6.4[ ). 
These estimants and their extrapolations are found in 



Table III . Note that we can no longer limit the parame- 
ter values by rib > as the condition (sb,s g ) = (1,0) 



treats the two flavours asymmetrically. In the case of the 
FPL model (rib = 1) it was found jlJJ that X was inde- 
pendent of n g . It is evident from our numerical data that 
this rig-independence in fact pertains to all rib € [0,2]. 
Final results for X as a function of rib have therefore 
been computed by averaging the available extrapolated 
scaling dimensions over n g . For rib = 1 the agreement 
with the result X w 1/8 found by Batchelor et al. |l4| ] 
is excellent. Furthermore we are able to conjecture the 
general formula, Eq. ( 7.15| ), for X as a function of the 
loop fugacities. 

When (sb,s g ) = (1,1) the parity of L must again 
be even, and we can make parabolic fits for the gap 



fo' 1] {L) - fo V '"'(L), as in Eq. (|10.5|) , without taking re- 
sort to the less accurate method of fitting for two central 
charges separately as above. The corresponding univer- 
sal amplitude is identified with the scaling dimension x\. 
The results, now for i max = 12, are shown in Table IV, 
and our values for the scaling dimension are once again in 
agreement with the analytical results, apart from n g = 2 
where logarithmic corrections are the most likely source 
of systematic errors [ p"l|| . 



f (o,o). 



rib 


n g 


X(3, 13) 


X(5, 13) 


X(7, 13) 


X(9, 13) 


Extrapolation 


Result 


Exact 


0.0 


0.0 


-0.05586 


-0.06109 


-0.06203 


-0.06232 


-0.06257 


-0.06269 (31) 


-0.06250 


0.0 


0.5 


-0.06080 


-0.06197 


-0.06220 


-0.06233 


-0.06253 




-0.06250 


0.0 


1.0 


-0.06043 


-0.06198 


-0.06221 


-0.06233 


-0.06250 




-0.06250 


0.0 


1.5 


-0.05869 


-0.06156 


-0.06215 


-0.06233 


-0.06259 




-0.06250 


0.0 


2.0 


-0.05804 


-0.06190 


-0.06297 


-0.06316 


-0.06324 




-0.06250 


0.5 


0.0 


0.04674 


0.04538 


0.04558 


0.04569 


0.04587 


0.04583 (16) 


0.04591 


0.5 


0.5 


0.04572 


0.04585 


0.04589 


0.04588 


0.04588 




0.04591 


0.5 


1.0 


0.04643 


0.04622 


0.04614 


0.04607 


0.04595 




0.04591 


0.5 


1.5 


0.04781 


0.04675 


0.04638 


0.04622 


0.04590 




0.04591 


0.5 


2.0 


0.04828 


0.04664 


0.04593 


0.04573 


0.04555 




0.04591 


1.0 


0.0 


0.11895 


0.12278 


0.12398 


0.12438 


0.12501 


0.12497 (8) 


0.12500 


1.0 


0.5 


0.12346 


0.12422 


0.12458 


0.12470 


0.12489 




0.12500 


1.0 


1.0 


0.12465 


0.12485 


0.12495 


0.12496 


0.12498 




0.12500 


1.0 


1.5 


0.12584 


0.12540 


0.12529 


0.12521 


0.12508 




0.12500 


1.0 


2.0 


0.12652 


0.12549 


0.12513 


0.12501 


0.12490 




0.12500 


1.5 


0.0 


0.17106 


0.18253 


0.18453 


0.18536 


0.18662 


0.18663 (25) 


0.18687 


1.5 


0.5 


0.18283 


0.18468 


0.18553 


0.18585 


0.18633 




0.18687 


1.5 


1.0 


0.18515 


0.18588 


0.18620 


0.18632 


0.18646 




0.18687 


1.5 


1.5 


0.18684 


0.18684 


0.18687 


0.18684 


0.18680 




0.18687 


1.5 


2.0 


0.18878 


0.18796 


0.18759 


0.18741 


0.18696 




0.18687 


2.0 


0.0 


0.2076 


0.2296 


0.2321 


0.2340 


0.2369 


0.2392 (27) 


0.2500 


2.0 


0.5 


0.2283 


0.2323 


0.2342 


0.2351 


0.2371 




0.2500 


2.0 


1.0 


0.2325 


0.2347 


0.2358 


0.2363 


0.2383 




0.2500 


2.0 


1.5 


0.2357 


0.2372 


0.2379 


0.2383 


0.2400 




0.2500 


2.0 


2.0 


0.2402 


0.2413 


0.2417 


0.2420 


0.2435 




0.2500 



TABLE III. Estimants X(Lo,L max ) for the scaling dimension of the twist operator along with their extrapolations to the 
infinite-system limit. For n b = 1 the value X = 1/8 was previously found to be independent of n s [14]. It is evident that this 
rig-independence holds for any value of rib, and in accordance herewith our final result is obtained by averaging the various 
extrapolants over n g . 
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Finally, the results for X2 as obtained from parabolic 
fits for the gap fo°\L) — /q (L) are shown in Table 
[v|. Again we have L max = 12. Just like in the case of X 
we find the extrapolated values of X2 to be independent 
of n s , and final results are obtained by averaging over 
this parameter. 



D. Entropy 

Apart from the various universal quantities, such as 
the central charge and the scaling dimensions, the trans- 
fer matrices also provide numerical values for the resid- 
ual entropy per vertex, s = /o(oo). In the limit n g — > 
of compact polymers this quantity is of interest to the 
protein folding community, due to the fact that native 
conformations of all globular proteins are compact S . 

Using our knowledge of the exact form of the finite-size 
corrections of order 1/L 2 , Eq. ( 10. 3| ), we have obtained 



very accurate extrapolations to the limit of an infinite 
system. 6 After subtracting the 1/L 2 correction a series 
of estimants s(L,L max ) may be obtained by fitting the 
residual size dependence to a pure 1/L 4 form. The re- 
maining L-dependence of these estimants turns out to be 
well accounted for by a further 1/L 4 fit, and in this way 
we arrive at a final value for s. The error bar on the final 
value can be estimated as its deviation from the most 
accurate extrapolant, s(X max — 2, L max ). 

The most accurate results are quite naturally found 
by employing this procedure on /q°'°^(L), and they are 



shown in Table VI. Results obtained by extrapolating 
the free energies for other sectors of the transfer matrix 
containing strings are consistent herewith but have er- 
ror bars that are roughly 10 times larger. If the fugacity 
of one of the strings equals two the error bars are even 
larger, which is to be anticipated from the fact that log- 
arithmic corrections to the scaling dimensions are larger 
than similar corrections to the central charge [pi]. 



n b 


n g 


zi(4,12) 


xi(6,12) 


a* (8, 12) 


Extrapolation 


Ref. [14] 


Exact 


0.0 


0.0 


-0.2433 


-0.2447 


-0.2470 


-0.2500 




-0.2500 


0.5 


0.0 


-0.1328 


-0.1295 


-0.1303 


-0.1313 




-0.1323 


0.5 


0.5 


-0.01713 


-0.01228 


-0.01217 


-0.01217 




-0.0131 


1.0 


0.0 


-0.0440 


-0.0423 


-0.0430 


-0.0439 


-0.0444 (1) 


-0.0446 


1.0 


0.5 


0.0737 


0.0763 


0.0764 


0.0765 


0.0750 (3) 


0.0761 


1.0 


1.0 


0.16608 


0.16646 


0.16657 


0.16663 


0.1667 (1) 


0.16667 


1.5 


0.0 


0.0267 


0.0271 


0.0264 


0.0255 




0.0260 


1.5 


0.5 


0.1466 


0.1472 


0.1472 


0.1472 




0.1483 


1.5 


1.0 


0.2411 


0.2395 


0.2395 


0.2394 


0.242 (2) 


0.2405 


1.5 


1.5 


0.3196 


0.3159 


0.3156 


0.3156 




0.3162 


2.0 


0.0 


0.0845 


0.0848 


0.0844 


0.0839 




0.1042 


2.0 


0.5 


0.2070 


0.2067 


0.2071 


0.2076 




0.2295 


2.0 


1.0 


0.3048 


0.3021 


0.3024 


0.3028 


0.307 (2) 


0.3250 


2.0 


1.5 


0.3882 


0.3841 


0.3842 


0.3843 




0.4044 


2.0 


2.0 


0.4640 


0.4618 


0.4635 


0.4657 




0.5000 



TABLE IV. Scaling dimension Xi, corresponding to one string of each flavour. 



In the special case of the equal-weighted six-vertex 
model, (rib,«g) = (1, 1), our value for s is in excellent 
agreement with the exact result due to Lieb p2[, 



a(l,l) = 




0.4315231 (10.7) 



and in the limit of two mutually excluding Hamiltonian 
walks, (n\y,n g ) = (0,0), we are able to conjecture the 
result 



s(0, 0) = - ln(2) ~ 0.3465735 ■ ■ • . (10.8) 



In fact, after having made this conjecture we discovered 
that the numerical values of /q 1 ' 1 '' (L), i.e., the free energy 
per site in the sector where we enforce one string of each 
flavour, are independent of L for 4 < L < 12, and equal 
to i ln(2) with full 16-digit machine precision. Since the 



The logarithmic corr ectio ns to the free energy implied by 
the A/" 7-1 term in Eq. (2.6) does not pertain to the cylindri- 
cal geometry implicit in our transfer matrix calculations. A 
similar remark applies to the surface term )l . 
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free energy per site in the thermodynamic limit is un- 
changed by the introduction of a string defect, this ob- 
servation lends credibility to the correctness of the above 
conjecture. 

Our result in the compact polymer limit merits spe- 
cial attention. Traditionally the entropy is quoted in 
terms of the so-called connective constant K — e^ 1 ' -*; 
see Eq. (2.6). Early approximations due to Flory |5^] 
and Huggins [B3 yielded respectively 



z — 1 

KFlory = — 1-104 



and 



K Huggins — ( z 1) I 1 



2V /2-i 3 



(10.9) 



(10.10) 



Here z = 4 is the coordination number of the square lat- 
tice. More recently, K ~ 1.472 was found from transfer 
matrix calculations [^6| and K — 1.475(15) by exhaustive 
computer enumeration of short-chain configurations [Q . 
Both these results are very close to the mean-field value 
kmf = - = 1.4715 ■ • • 0, and it is tempting to conclude 
that conformations of compact polymers are in fact de- 
scribed by mean-field theory |4j. However, our result 



k = 1.472801(10) 
demonstrates that this is not the case. 



(10.11) 



XI. DISCUSSION 

From the construction of the effective field theory of 
the FPL 2 model some rather general conclusions regard- 
ing the scaling of compact polymers, and the relation 
between loop models and conformal field theory can be 
drawn. It also provides new insights into the three-state 
Potts antiferromagnet and the dimer loop model, which 
are identified with specific points in the phase diagram 
of the FPL 2 model. Wc conclude the paper with a dis- 
cussion of these topics. 



A. Compact polymers 

One of the main motivations for studying fully packed 
loop models is provided by compact polymers, their scal- 
ing properties in particular. Just like polymers in the 
dilute and dense phase, compact polymers form a crit- 
ical geometrical system characterised by conformational 
exponents 7 and v. The exponent 7 relates the num- 
ber of conformations of the polymer to the number of 
monomers; see Sec. [il] for details. The other conforma- 
tional exponent (y) relates the linear size of the polymer 
to the number of monomers. For compact structures it 
has the trivial value 1/2 since these polymers are space 
filling. 



Tib 


n g 


Z2(4, 12) 


z 2 (6,12) 


1-2(8,12) 


Extrapolation 


Result 


Exact 


0.0 


0.0 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 (0) 


0.0000 


0.0 


0.5 


0.0000 


0.0000 


0.0000 


0.0000 




0.0000 


0.0 


1.0 


0.0000 


0.0000 


0.0000 


0.0000 




0.0000 


0.0 


1.5 


0.0000 


0.0000 


0.0000 


0.0000 




0.0000 


0.0 


2.0 


0.0000 


0.0000 


0.0000 


0.0000 




0.0000 


0.5 


0.0 


0.1279 


0.1355 


0.1372 


0.1389 


0.1386 (2) 


0.1386 


0.5 


0.5 


0.1365 


0.1371 


0.1378 


0.1387 




0.1386 


0.5 


1.0 


0.1377 


0.1374 


0.1379 


0.1385 




0.1386 


0.5 


1.5 


0.1383 


0.1375 


0.1379 


0.1384 




0.1386 


0.5 


2.0 


0.1392 


0.1376 


0.1379 


0.1383 




0.1386 


1.0 


0.0 


0.2333 


0.2447 


0.2472 


0.2504 


0.2495 (5) 


0.2500 


1.0 


0.5 


0.2488 


0.2477 


0.2484 


0.2493 




0.2500 


1.0 


1.0 


0.2514 


0.2487 


0.2490 


0.2494 




0.2500 


1.0 


1.5 


0.2538 


0.2497 


0.2495 


0.2492 




0.2500 


1.0 


2.0 


0.2573 


0.2512 


0.2504 


0.2494 




0.2500 


1.5 


0.0 


0.3197 


0.3377 


0.3416 


0.3466 


0.3487 (26) 


0.3506 


1.5 


0.5 


0.3429 


0.3425 


0.3443 


0.3466 




0.3506 


1.5 


1.0 


0.3486 


0.3457 


0.3466 


0.3478 




0.3506 


1.5 


1.5 


0.3548 


0.3497 


0.3496 


0.3494 




0.3506 


1.5 


2.0 


0.3636 


0.3561 


0.3547 


0.3529 




0.3506 


2.0 


0.0 


0.3920 


0.4202 


0.4268 


0.4353 


0.446 (12) 


0.5000 


2.0 


0.5 


0.4244 


0.4277 


0.4323 


0.4382 




0.5000 


2.0 


1.0 


0.4346 


0.4348 


0.4382 


0.4426 




0.5000 


2.0 


1.5 


0.4468 


0.4452 


0.4474 


0.4502 




0.5000 


2.0 


2.0 


0.4640 


0.4618 


0.4635 


0.4657 




0.5000 


TABLE V. Scaling dimension 12, corresponding to two black strings. 
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Prior to our work, exact results have been obtained 
for compact polymers on the Manhattan |35| and the 
honeycomb |20| lattice, and the mean- field value 7=1 
was found in both cases. This value of 7 indicates that 
the two ends of the compact polymer are independent 
at large distances. This follows from the scaling rela- 
tion x\ = 1 — 7 = 0, where the one-string dimension 
xi describes the probability G\{r) ~ r ~ 2xi that the two 
chain ends are separated by a distance r. In this re- 
gard the scaling of compact polymers on the Manhattan 
and the honeycomb lattices is equivalent to that of ideal 
chains. Ideal chain configurations are described by sim- 
ple random walks for which each step is independent of 
the previous one. 

Here we have calculated the exact conformational ex- 
ponent 7 = 117/112 for compact polymers on the square 
lattice. The fact that 7 > 1 is tantamount to an effective 
repulsion between the ends of the chain, indicating non- 
ideal behaviour. Indeed, the fact that the connective 



constant k in Eq. ( 10.11 ) is larger than its mean- field 
value indicates that the origin of this repulsion is en- 
tropic. Earlier numerical studies of this problem utilising 
direct enumerations of chain conformations have failed to 
see any deviation from the ideal chain result 7mf = 
we can attribute this to the fact that the actual differ- 
ence is indeed very small (7 — 7mf = 5/112) and below 
the numerical accuracy previously achieved. The same 



comment can be made for the connective constant. 

Another interesting aspect of compact polymers is that 
their scaling properties are lattice dependent. This is 
in contrast to the dilute and dense case which are de- 
scribed by conformational exponents that do not depend 
on the lattice type (e.g., honeycomb versus square). As 
remarked earlier this "lack of universality" is due to a 
kind of geometrical frustration that arises from the fully 
packing constraint imposed on the loop models which are 
employed in studies of compact polymers. 

Finally, the field theory solution of the FPL 2 model 
uncovered a property of compact polymers that, to our 
knowledge, was not previously anticipated. The fact that 
there is a whole line of critical points in this loop model in 
the Hamiltonian walk limit (ny, —> 0) indicates a contin- 
uum of universality classes described by compact poly- 
mers on the square lattice. In particular the exponent 
7 can be changed continuously by adjusting the fugac- 
ity of the loops uncovered by the polymer. The loop 
weight of the uncovered (grey) loops can be thought of 
as an effective interaction amongst the monomers, albeit 
a non-local one. A similar effect of interactions on di- 
rected self-avoiding walks was discovered by Cardy |3(i| l 
from a field theoretical analysis of the problem. The exis- 
tance of a continuously varying 7 in this case was recently 
challanged by numerical results [plj . 



nb 


n g 


/r } (4) 


/r»(6) 


/r>( 8 ) 


/r>(io) 


AT (12) 


zr>(i4) 


s 


0.0 


0.0 


0.17328680 


0.28881133 


0.31784496 


0.32923359 


0.33490107 


0.33815371 


0.346575 (14) 


0.5 


0.0 


0.26740000 


0.33317928 


0.35057672 


0.35745438 


0.36088114 


0.36284872 


0.367950 (9) 


0.5 


0.5 


0.35063553 


0.37668215 


0.38371283 


0.38639599 


0.38769210 


0.38842126 


0.390258 (3) 


1.0 


0.0 


0.32923947 


0.36764369 


0.37752555 


0.38137032 


0.38327066 


0.38435762 


0.387166 (7) 


1.0 


0.5 


0.40772622 


0.41103439 


0.41126990 


0.41111188 


0.41095017 


0.41082815 


0.410405 (2) 


1.0 


1.0 


0.46298939 


0.44576535 


0.43960110 


0.43671524 


0.43513763 


0.43418273 


0.4315233 (4) 


1.5 


0.0 


0.37601935 


0.39599984 


0.40063320 


0.40233073 


0.40314475 


0.40360330 


0.404771 (5) 


1.5 


0.5 


0.45180855 


0.43964968 


0.43509788 


0.43291844 


0.43171625 


0.43098591 


0.4289459 (10) 


1.5 


1.0 


0.50624745 


0.47501911 


0.46431698 


0.45948057 


0.45688890 


0.45533728 


0.4510742 (17) 


1.5 


1.5 


0.54930614 


0.50513652 


0.49006459 


0.48331974 


0.47972832 


0.47758588 


0.471726 (2) 


2.0 


0.0 


0.41389271 


0.42018005 


0.42097629 


0.42111147 


0.42113891 


0.42114428 


0.421145 (6) 


2.0 


0.5 


0.48795109 


0.46429984 


0.45622604 


0.45257417 


0.45061901 


0.44945033 


0.4462607 (10) 


2.0 


1.0 


0.54202495 


0.50046092 


0.48641918 


0.48016010 


0.47683419 


0.47485271 


0.4694505 (18) 


2.0 


1.5 


0.58515036 


0.53158535 


0.51333087 


0.50520022 


0.50088581 


0.49831761 


0.491323 (3) 


2.0 


2.0 


0.62122666 


0.55918707 


0.53795845 


0.52850379 


0.52348906 


0.52050483 


0.5123870 (19) 



TABLE VI. Residual entropy s, obtained by extrapolating /q (L) to the infinite-system limit. 



B. Relation to other models 



1. Dimer loop model 



The FPL 2 model is a loop model which exhibits a two- 
dimensional manifold of fixed points in its phase diagram. 
Certain points in the critical region map to previously 
studied lattice models and here we comment on the rel- 
evance of our results for these models. 



The dimer loop model studied by Raghavan et al. p| 
is the n\y — 2, Hg = 1 FPL 2 model; see Fig. |l|. The dimer 
loop model is defined by placing black and white dimers 
on the square lattice so that every vertex of the lattice is 



25 



covered by exactly one black and one white dimer. Every 
such configuration is given equal weight. The mapping to 
the FPL 2 model is achieved by identifying the bonds cov- 
ered by dimers as making up the black loops, whilst the 
uncovered bonds form the grey loops. The original mo- 
tivation for studying this dimer problem is that it leads 
to a height model with a two-component height; cfr. the 
traditional dimer model which is described by a single 
component height. 

Performing Monte Carlo simulations of the dimer loop 
model Raghavan et al. reached the conclusion that one 
of the two height components is rough whilst the other 
one is "anomalously smooth", i.e., its structure function 
decays at small wave- vectors q slower then 1/q 2 ; a 1 /q 2 
dependence is to be expected in a Gaussian field theory. 

In light of our results we would conclude that the dimer 
loop model is criti cal with a central charge c = 2. This 
follows from Eq. (OA) for rib = 2 and n g = 1. The 
two components of the height found by Raghavan et 
al. should therefore both be rough, each contributing one 
to the central charge (c = 1 + 1). Furthermore, we be- 
lieve that the observed anomalous behaviour of one of 
the heights can be attributed to the fact that this model 
is exactly at the boundary of the critical region of the 
FPL 2 model. We observe a similar effect in our numeri- 
cal transfer matrix results which show largest deviations 
from the proposed exact formulae for loop fugacities at 
the critical-region boundary. The culprit might be loga- 
rithmic corrections due to the presence of marginal oper- 
ators. To check this hypothesis and reconcile it with the 
fact that no such effects are seen in Monte Carlo simu- 
lations of the four-colouring model p^] (rib = n g = 2), 
which is also at the boundary of the critical region, sim- 
ulations of the dimer-loop model for larger system sizes 
would be welcome. 



2. Three-state Potts antiferromagnet 

The critical ground state of three-state Potts antiferro- 
magnet maps to the equal- weighted six- vertex model p"6| ] 
which is the rib = n s = 1 point in the critical region of 
the FPL 2 model; see Fig. [j]. Along the line rib = n g the 
colouring representation of the FPL 2 model has the ad- 
ditional symmetry with r espect to cyclic permutations of 
the four colours; see Sec. V C 2 . This explains the origins 
of the Z4 symmetry found by Salcur for the three-state 
Potts antiferromagnet [p6|. 



3. Folding model 

The folding model of the square-diagonal lattice re- 
cently investigated by Di Francesco j57| maps onto a 
constrained version of the (rib, fig) = (2,2) FPL 2 model. 
The constraint consists in allowing only the vertices 1, 
3, 5 and 6 of Fig. || for sites on the even sublattice, and 
similarly vertices 2, 4, 5 and 6 on the odd sublattice. 

We have modified our transfer matrices to take this 
constraint into account. Our result for the folding en- 
tropy, s — 0.4604(4), is in complete agreement with 
Rcf. ^7|. 7 Interestingly enough the finite-size scaling of 
the gaps in the eigenvalue spectrum seems to indicate 
that the model is not critical for general values of the 
loop fugacities. From the field theory of the FPL 2 model 
we should be able to to understand why the constraint 
imposed by the folding model leads to a relevant pertur- 
bation which takes the system away from criticality. This 
we leave as an interesting open question. Incidentally, 
the situation is very reminiscent of the reformulation of 
the Q-state Potts model in terms of a staggered vertex 
model. Only at the critical point are the vertex weights 
on the even and odd sublattices identical, thus allowing 
for an exact solution of the model fl32|] . 



C. Conformal field theory 

The Liouville field theory proposed for the effective 
theory of the FPL 2 model in the critical region is confor- 
mally invariant. Each point in the critical phase is char- 
acterised by the central charge and the scaling dimen- 
sions of primary fields, which are associated with electric 
and magnetic charges in the Coulomb gas. For generic 
values of the loop fugacities the background charge e is 
not commensurate with the electric charges that make 
up the lattice TZ* . This implies that amongst the electric 
operators there will be many (an infinite number, in fact) 
that have negative dimensions, signaling the non- unitary 
nature of the conformal field theory. Non- unitary CFT's 
appear in many other critical geometrical models, critical 
percolation being the best known example. 

Liouville field theory provides the Euclidean action for 
the Coulomb gas description of conformal field theories 
proposed by Dotsenko and Fateev Q. As such it con- 
tains the so-called screening charges which are the ver- 
tex operators that make up the Liouville potential. In 
the original formulation these charges were introduced 
on formal grounds so as to ensure the existence of non- 
vanishing four-point correlation functions in the theory. 



7 Our normalisation is "per vertex" whilst that of Di 
Francesco is "per triangle". Accordingly we find twice his 
result. 
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In order for these vertex operators not to disrupt the con- 
formal symmetry of the modified Gaussian model (the 
modification is the addition of the boundary term to 
the gradient-square action) they are necessarily marginal, 
i.e., their scaling dimension is 2. 

Here we have found a physical interpretation of the 
screening charges. Their role in loop models is to ensure 
that the number of large loops from scale to scale stays 
of order one; this translates into the statement that the 
loop fugacities do not flow under the action of the renor- 
malisation group. 

The fact that we have a concrete physical interpreta- 
tion of the screening charges directly leads to the calcula- 
tion of the elastic constants in the Liouville field theory. 
In the traditional Coulomb gas approach these coupling 
constants are calculated by comparing to formulae de- 
rived from an exact solution of the model. Once these 
constants are known marginal vertex operators that play 
the role of screening charges can be written down. Our 
construction basically reverses this procedure, and by do- 
ing so makes no reference to an exact solution. 

Finally, we end with a speculatory note concerning the 
prospects of solving the FPL 2 model via Bethe Ansatz. 
Namely, all loop models to date have been solved by this 
method after mapping them to a vertex model, follow- 
ing a procedure analogous to the one outlined in Sec. IV. 
This does not seem to work for the FPL 2 model, at least 
not along the n g = 1 line Why this is so is an 

interesting open question. 

One possibility is that the full FPL 2 model needs to 
be considered as opposed to the FPL model studied by 
Batchelor et al. for which n g — 1 is fixed. A more intrigu- 
ing possibility is that a Bethe Ansatz solution might be 
hindered (or made more difficult) by the non-trivial elas- 
ticity displayed by the FPL 2 model in its interface rep- 
resentation. This statement we base solely on the obser- 
vation that all previously solved loop models are simple 
as interface models in the sense that the height fluctua- 
tions are described by a single elastic constant. For the 
FPL 2 model, as described in Sec. 0, the stiffness tensor 
consists of three independent components. Whether in- 
deed the interface representation of the loop model has 
any bearing on its Bethe Ansatz solvability remains to 
be seen. 
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APPENDIX A: DIMENSIONS OF ELECTRIC 
AND MAGNETIC OPERATORS 

We calculate the scaling dimensions of electric and 
magnetic operators in the Coulomb gas theory described 
by the action 

S cg = IJ rf 2 x 9a (dH a ) 2 + JL />x (Eo • H)K , 

(Al) 

where TZ is the scalar curvature. We are interested in 
the situation when the height field is defined on a flat 
surface, in which case TZ is zero everywhere except at the 
boundaries. 



1. Electric charges 

The scaling dimension x(E), of the electric-type oper- 
ator exp(iE • H(x)), follows from the two-point function 



iE-H(x) -t(E-2E )-H(y) 



(A2) 



where the expectation value is with respect to the mea- 
sure defined by the action Sec ■ The extra electric charge 
2Eo appears due to the charged boundary conditions en- 
forced by the curvature term in the Coulomb gas action, 
Eq. @. 

We break up the calculation into two parts. First we 



calculate the two-point function, Eq. (A2), in the absence 
of the background charge (Eo = 0). We make use of the 
property of Gaussian integrals, 



e iE.H(x) e -iE-H(y) 



(A3) 



exp(--(i; Q ) 2 ((i^(x)-i^(y)) 2 ) 



and of the known propagator for the massless scalar field 
in two dimensions (where we have dropped the regulators 
at large and small distances), 

((H a (x)-H<*(y)) 2 ) = — ln|x-y| . (A4) 

Combining the above two equations and comparing the 
result to Eq. ([A3), we find 



2^ 0) (E) = 



1 



2irg a 



(A5) 



the superscript (0) is there to remind us that this formula 
is valid only for Eq = 0. 
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This result for the two-point function can be rewritten 



as 



where 



e ^H( x ) e - iE .H (y )\ =exp[f (0) (X)y)L (Ag) 



4° ; (x,y) = - 



2Trg a 



(£ Q ) 2 ln|x-y| 



(A7) 



is the energy for two (vector) electric charges interacting 
via the two-dimensional Coulomb force; in this language 
«Se is the energy of the electrostatic field set up by the 
electric charges ±E, expressed in terms of the electro- 
static potential h. This seemingly trivial rewriting makes 
the calculation of x(E), the electric dimension in the pres- 
ence of a background charge, physically transparent. 

To properly take into account the curvature term we 
define the height field over a disc of radius R, instead 
of the infinite plane, keeping in mind that at the end 
of the calculation we need to take the limit R — > oo. 
In the case of the disc 1Z = 8ttS(R), and the curva- 
ture term introduces a charge — 2Eo at the disc bound- 
ary. Therefore, the vacuum of the modified Coulomb 
gas must contain a floating charge +2Eo in the disc in- 
terior, and the electrostatic energy of this charged vac- 
uum is £^ o (0,R) = ~AEl a \n(R)/2ng a . Now, to find 
the scaling dimension of a vertex operator of charge E, 
we imagine placing charges +E and — E at points x and 
y in the disc interior, and we calculate the total electro- 
static energy with respect to the charged vacuum. The 
floating charge being positive will coalesce with the neg- 
ative charge — E. Using Coulombs law, Eq. (A7), we 
then calculate the interaction energy of charges +E at 
x, E + 2Eo at y, and — 2Eo at R, keeping in mind 
R > |x - y|. The final result 



£e(x, y) 



1 



2irg a 



E a (E a -2E 0a )ln\x-y\ (A8) 



is obtained after the energy of the charged vacuum is 
subtracted. Now it is a simple matter to read off the 
scaling dimension as the negative coefficient in front of 
the logarithm, 



2x(E) 



1 



■E a {E a — 2Eo a ) 



(A9) 



This result can be derived in a more rigorous fashion by 
constructing the stress-energy tensor for the field theory 
Scg and calculating its operator product with the vertex 
operator exp(iE • H) |3q]. 



Z >M (x,y)/Z 5 



-2x(M) 



(A10) 



■Z>m(x, y) is the sum (path integral) over height config- 
urations where a vortex and an antivortex, of topological 
charge ±M, are placed at positions x and y of the basal 
plane, whilst Z > is the unconstrained sum: 

Z > = JvH cxp f-i J rf 2 x g a (dH a ) 2 ^j . (All) 

Here we have dropped the curvature term since it does 
not affect correlation functions of magnetic operators. 

We can use the electrostatic analogy once again. 
Namely, we consider the interaction energy between two 
topological defects, £m(x, y). Since Z > is a Gaussian 
path integral, it follows that 



Z>M(x,y)/Z> = exp[£ M (x,y)] 



where 



£ivr(x,y) 



2tt 



(M Q ) 2 ln|x-y| 



(A12) 



(A13) 



The above interaction energy is calculated as the Gaus- 
sian action of the the classical configuration of the height 
field, h c . h c solves the classical equations of motion 
(Laplace's equation) with boundary conditions dictated 
by the presence of topological defects at x and y 63] . The 
scaling dimension of a magnetic-type operat or is then the 
coefficient in front of the logarithm in Eq. ( A13), 



2x(M) 



9a_ 

2ir 



{M a f 



(A14) 



APPENDIX B: ENUMERATION OF THE 
CONNECTIVITIES 

The implementation of the transfer matrix (TM) for 
the FPL 2 model on a cylinder of width L and length 
M requires an enumeration of the possible connectivity 
states of the L points on the dangling edges of row M. 
Each of these L points can either 

1. be connected by Gm to one of the dangling edges 
of row through a string of flavour i = b, g, or 

2. be connected by Gm to one and only one other point 
in row M through a loop segment of flavour i — b, g. 

A suitable representation of this information is furnished 
by a double state vector 



2. Magnetic charge 

To calculate the magnetic dimension x(M) we consider 
the ratio of partition functions, 



6 l t 2 6 3 
l l l 2 l 3 ' 



(Bl) 



which we shall refer to as the index representation. The 
indices i\ (k = 1,2, .... L) are defined as follows: 
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1- if = if is a (non- unique) positive integer if and 
only if points k and I are interconnected through a 
black string. 

2. if = if and only if point k touches a grey string 
or loop segment. 



the relevant sector. The practical implementation of the 
TMs relies on the mapping from the index to the number 
representation and its inverse. 

We shall now consider, one by one, the various sectors 
of the TM. 



3. i\ 



T if and only if point k is connected to a 



dangling edge of row through a black string. 

A similar definition is true for the indices if provided 
that one reads "grey" instead of "black" and vice versa. 
Two index representations are said to be identical if they 
are so up to the arbitrariness of the choice of positive 
integers. Also note that if if ^ we have if = and 
conversely. 

A restriction on those indices that take positive val- 
ues follows from the fact that loops of the same flavour 
are not allowed to intersect. Namely, if j < k < I < m 



the equalities i° = if and if 



cannot both be true. 



So in addition to being pairwise these connectivities are 
also well-nested p5| . The same is true for the grey in- 
dices, whereas there are no such restrictions when both 
flavours are involved. Indeed, connectivity states with 
i}j = if and if = if n are explicitly allowed by the last two 
vertices shown in Fig. |^. 

In practice we are only interested in the first few eigen- 
values of TMs having a definite number of strings of 
each flavour. The relevant sectors of the TM are de- 
noted T( Sb ' Sg ), where s, is the number ofstrings of flavour 
i = b, g. The fully packing constraint means that we can 
only examine system sizes L that have the same parity 
as Sb + s g- The various sectors have different physical 
interpretations and each requires a different enumeration 
of the connectivity states. Since the two flavours enter at 
an equal footing in the partition function, Eq. (9.1), wc 
only need consider Sb > s g . The T^ 0,0 ) sector contains 
information about the free energy and the energy-like 
correlation length. The geometrical scaling dimensions 
Xi and X2 can be obtained from the T^ 1,1 ) and the T^ 2,0 ) 
sectors respectively. Finally the sector T^ 1,0 ) gives the 
scaling dimension of the twist-like operator. 

Whilst the index representation contains all informa- 
tion necessary for determining the value of a given entry 
in the TM it is obviously not suitable for labeling the 
entries. We therefore need another representation, the 
so-called number representation, in which the connectiv- 
ities are labeled by the integers 1,2, 



.,c! Sb,5s ', where 



C 



1. T 



(0,0) 



sector 



When no strings are present all the L dangling edges 
of row M are pairwise connected with either a black or 
a grey loop segment. In particular L must be even. For 
any particular connectivity we can then decompose L as 
L = 2pb + 2p g , where pi > is the number of pairs of 
dangling edges covered by a flavour i loop segment. Since 
loops of different flavours are allowed to cross (see Fig. |^) 
the total number of connectivities is 



C 



(0,0) _ 



L=2p b +2p g 



L 
2p, 



(B2) 



where c p is the number of pairwise well-nested connec- 
tivities of 2p points. The c p 's were first considered in 
the context of the Potts model |45|l , but were also found 
to play a central role in the TM formulation of the O(n) 
model [[l7| . We shall now briefly recall how they are eval- 
uated. 

Consider a well-nested pairwise connectivity of 2p 
points given by the index representation (iii2 . . . i2p)- A 
recursion relation follows from observing that i\ = i2k 
for precisely one integer k > 1. According to the well- 
nestedness criterion the sub-sequences (i2*3 ■ • • i2fc-i) and 
(i2fc+ii2fe+2 • • • iip) are both well-nested, and indices oc- 
curring in one of them do not occur in the other. Hence 
for p > 1 



fc=i 



(B3) 



and Co = 1. By means of the generating function 
P{ x ) = J2"^Lo c p xP it ^ s rea dily shown Q that 



(2p)\ 



P \(p + l)V 



(B4) 



and that asymptoticall y c v ~ 4 P . 

Using Eqs. (B2) and (B4) we can now compute explicit 



values for the C£'°^ ■ These are shown for 2 < L < 16 in 



is the number of different connectivity states in Table VII. 



L 


A L 


c (o,o) 




(^(2,0) 


L 


C£.0) 


2 


16 


2 


2 


1 


1 


1 


4 


256 


10 


24 


12 


3 


6 


6 


4,096 


70 


300 


150 


5 


50 


8 


65,536 


588 


3,920 


1,960 


7 


490 



29 



10 


1,048,576 


5,544 


52,920 


26,460 


9 


5,292 


12 


16,777,216 


56,628 


731,808 


365,904 


11 


60,984 


14 


268,435,456 


613,470 


10,306,296 


5,153,148 


13 


736,164 


16 


4,294,967,296 


6,952,660 


147,232,800 


73,616,400 


15 


9,202,050 



TABLE VII. The number Cj^ h ' Sg of FPL 2 connectivity states for L dangling edges accommodating Sj strings of flavour 
i = b, g. Only values of L with the same parity as Sb + s g are shown. When more than one string of any flavour is present 
further restrictions than the well-nestedness criterion apply, as described in the text. Accordingly the number C^ 2 ' ' is merely 
a useful upper limit on the true C} 2 ' '. The efficiency of writing the TMs in the connectivity basis can be appreciated by 
comparing C ( °' 0) to 4 L , the latter being the dimensions of the TM written in the conventional colour basis, where every 
dangling end is labeled independently by A, B, C or D. 



For obvious reasons we shall call the function 



p{i\h ■ ■ -h p ) = k 



(2.5) 



defined by i\ = 12k the cut function of the index repre- 
sentation (iiia ■ ■ ■ *2p)- A complete ordering of the well- 
nested sequences is now induced by applying the cut func- 
tion first to the whole sequence, then recursively to its 
right and finally to its left part |4^,|1J . Accordingly, the 
mapping from the index to the number representation for 
a well-nested one-flavour connectivity is accomplished by 



a{ii%2 ■ ■ ■ «2 P ) 



1 if p < 1 

Yll=l Cl-lCp-l + a(%2 ■ ■ ■ «2fc-l) 

+ [cr(i2fc+i • • • i-2 P ) - l]cfc_i otherwise, 

(2.6) 



where the c p are given by Eq. (B4). 

To give a complete specific atio n of the connectivity of 
any one flavour in the state (|Bl| ) we need to keep track 
of the positions of those indices that are zero. For a fixed 
number of z zero indices this is accomplished by the lex- 
icographic ordering 



1 if L-- 

= { ^(*2«3 • 

(V) - 



: 1 or z = L 

..tzOifii^O (2.7) 
■ipfaia ■■■ih) if h = 0, 



assigning the lowest value to the sequence with all the 
zeros accumulated to the right. 

The number representation of the two-flavour state 



(Bl) is now obtained by first ordering according to the 



number of indices vj. being zero, then lexicographically 
ordering the positions of these zero indices, and finally 
using the ordering ( |2.6| ), first on the well- nested subse- 
quence of non-zero black indices and then on the cor- 
responding grey subsequence. More precisely, the map- 
ping from the index to the number representation in the 
(sb, s g ) = (0, 0) sector is given by 



2-fc 



(2- 



+ [^(i b ) - l]c Pb c p& + [a(i b ) - l]c Pb + 



where i = (i^i^ ■ ■ ■ denotes the sequence of black in- 
dices and i b the subsequence of the pb pairs of non-zero 
indices (and, of course, similarly for the grey flavour). 

The inversion of Eq. ( |2.9| ) , so as to furnish a mapping 
from the number to the index representation, is straight- 
forward if we know how to invert the functions a and ip. 
Details on this have already been given in Ref. M ■ 



B. T (1 ' 0) sector 



In the case of one black string spanning the length of 
the cylinder the number of dangling edges in row M can 
be written as L = 2p^ + 2p g + 1, where the Pi have the 
same meaning as above. In particular L must be odd. 

The presence of one string of cither flavour does not im- 
pose any additional restrictions on the connectivity states 
of the subsequence of positive indices of that flavour. In- 
deed, if the position of the string is given by i b = — 1 



the non-zero subsequence of (i b +1 



_ 1 ) is still 



well-nested, and the arguments given above apply. The 
number of connectivity states is therefore found by mul- 
tiplying the L possible positions of the string by the num- 
ber of (sb,s g ) = (0,0) states of the remaining L—l points 



c (i,o) = L q(o,o) 



(2.9) 



Explicit values are shown in Table VII 



Similarly the mapping from the index to the number 
representation is found by first ordering after the position 
r of the string, and then after the value of <f>(°>°> taken of 
the remaining indices 



(1,0) 



,-b„-b 



(r - 1)C 



(0,0) 
L-l 



(2.10) 
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C. T 



sector 



When one string of each flavour is present L = 2pb + 
2p g + 2 must be even, and again it suffices to augment the 
considerations from the T' 0,0 ) case by some book-keeping 
as to the positions of the two strings. Explicit values of 



C 



L(L - 



(2.11) 



are shown in Table VII 



Letting r, denote the position of the string of flavour 
i = b, g we find that 



(2,0) 



.0(0,0) 



= [(rb-l)(L-l) 



,b -b 

■ Vb-iVb+i • ■ 

■ V b -iV b +i ■ ■ 



AS 

-iV+i • 



(2.12) 



is the desired mapping from the index to the number 
representation. 

A possible configuration of the system for (sh,s g ) = 
(1, 1) is illustrated in Fig. |8|, where the index representa- 
tion of the connectivity state for each completed row is 
shown to the right of the figure. 



D. T 



(2,0) 



sector 



Considering now the case of two black strings, it ap- 
pears that the number of connectivity states for L even 
is given by 



C 



(2,0) 



C 



(0,0) 
L-2 ' 



where we have simply divided Eq. ( |2.11 ) by 2 to take 
into account the indistinguishability of two strings of the 
same flavour. This is however not quite true, since for 



L > 4 the number (2.13) includes certain disallowed basis 
states. For L = 4 these are 



-11-11 




and 



1-11-1 




(2.14) 



The reason why these states are not valid is that, by 
definition of the allowed vertices (see Fig. , black loop 
segments cannot cross a black string. In general, there- 
fore, any configuration where the positions of two equal, 
positive black indices are separated by exactly one black 
string is not a valid one, even though the positive in- 
dices of each flavour satisfy the well-nestedness criterion. 
Accordingly, the true cjf' ^ is less than the C^' ' of 



of L in order to determine the scaling dimensions x\ and 
X2 with the same numerical precision. Instead we found 



it efficient to construct all the C 



(2,0) 



basis states, list the 



number representations of those that are disallowed, and 
force the corresponding entries of T^ 2,0 ) to zero. 

With this proviso the mapping from the index to the 
number representation is 



0(2,0) ( \ 



l.O) 



if.. 



'ri-l'ri + l ' 



'rs-l'Va + l • ' 



(2.15) 



where r± and r-i are the positions of the two black strings, 
and tp(i h + 1) means that we should lexicographically or- 
der the positions of the black indices that are —1. 
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